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Dieses Buch stimmt in der Gliederung des Stoffes 
und in der Bezeichnungsweise überein mit dem uni-text 


Wolfgang Böhm und Günther Gose 
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Numerischen Mathematik 

1977. VIII, 152 Seiten 


„In dem Buch werden die Grundideen numerischer Lösungsmethoden 
dargestellt. Es ist für Studenten der Mathematik und der Informatik 
und für Interessenten naturwissenschaftlicher Disziplinen geschrieben. 
Nach einer Einführung von Grundbegriffen (Algorithmen, Matrizen) 
werden im zweiten Teil numerische Fragen der linearen Algebra be- 
handelt, darunter der Gaußsche Algorithmus mit Pivotsuche, der 
konzentrierte Gauß-Algorithmus, die Methode von Cholesky zur 
Erhaltung der Symmetrie, die Relaxationsmethode, die angenäherte 
Lösung über- oder unterbestimmter linearer Gleichungssysteme und 
Jas Simplexverfahren der linearen Optimierung. Iterative Verfahren 
wie die klassische Vektoriteration, der Rutishauser-Algorithmus, 
Methoden zur Konvergenzbeschleunigung und Nullstellenbestimmung 
bilden den dritten Teil. Es folgen die Interpolation durch Polynome 
und Splinefunktionen. 

Im abschließenden fünften Teil werden die numerische Differentiation 
und Integration einschließlich Fehlerabschätzungen, die Grenzwert- 
bestimmung durch Extrapolation sowie Ein- und einfache Mehrschritt- 
verfahren mit Schrittweitensteuerung zur Lösung von Differential- 
gleichungen behandelt, dabei wird auf das Runge-Kutta-Verfahren 
besonders eingegangen und ein Vergleich der Vorteile von Ein- und 
Mehrschrittverfahren durchgeführt. 

Die Darstellung strebt mehr nach Klarheit der Grundideen als nach 
Vollständigkeit oder weitestgehender Allgemeinheit. Die einzelnen 
Verfahren werden durch Algorithmen ergänzt, die die Formulierung 
in der Programmiersprache ALGOL vorbereiten. Das ausgezeichnete 
Lehrbuch ist besonders auch zum Selbststudium geeignet.” 
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Vorwort 


In den letzten Jahren war in der Herstellung immer leistungsfähigerer programmierbarer 
Taschenrechner eine rasante Entwicklung zu beobachten. Um ihre Möglichkeiten und 
Kapazitäten optimal auszuschöpfen, sollten auch für diese Kleinrechner Programmbiblio- 
theken zur Verfügung stehen. 

Der vorliegende Band enthält eine Sammlung nützlicher Programme der numerischen 
Mathematik für den programmierbaren Taschenrechner TEXAS INSTRUMENTS TI 59. 
Zugrundegelegt wurde das im gleichen Verlag erschienene Buch 


„Einführung in die Methoden der Numerischen Mathematik’' 


von Wolfgang Böhm und Günther Gose (Vieweg, Braunschweig 1977), aus dem Gliederung 
und Bezeichnungsweise übernommen wurden, um die Anwendung und das Arbeiten mit den 
Programmen zu erleichtern. Hier findet der interessierte Leser neben der theoretischen 
Herleitung auch die flußdiagrammähnlichen Algorithmen, nach denen die Programme 
erstellt wurden. 

Einen kurzen Überblick über die Handhabung des Rechners liefert das Kapitel O „Einführung”. 
Dem im Umgang mit dem TI 59 ungeübten Leser sei zunächst ein intensives Studium der 

zum Rechner gehörenden Bedienungsanleitung ‚Individuelles Programmieren”, insbesondere 
der Abschnitte I, Il und VII empfohlen. 

Mein ganz besonderer Dank gilt Herrn Prof. Dr. Wolfgang Böhm, dessen Vorlesungen ich die 
Freude an der numerischen Mathematik verdanke. Ohne seine Anregungen und aufmunternden 
Ratschläge wäre dieses Buch nicht entstanden. 

Nicht zuletzt danke ich der Firma TEXAS INSTRUMENTS für die freundliche Unterstützung und 
dem Vieweg Verlag für die problemlose Zusammenarbeit. 


Jürgen Kahmann 


Wolfenbüttel, im Frühjahr 1980 
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O Einführung 


Dieses Buch enthält 44 Programme der numerischen Mathematik für den programmierbaren 
Taschenrechner TEXAS INSTRUMENTS TI 59. Die ausführlichen Prograrmmbeschreibungen 
umfassen jeweils 


— eine kurze Erläuterung des programmierten Verfahrens 

— eine Tabelle der bei der Benutzung des Programms durchzuführenden Tastenfolgen 

— eine Übersicht über die Registerinhalte 

— ein vom TI 59 durchgerechnetes Beispiel und 

— einen vollständigen Programmausdruck (sowie bei den Programmen 4.3 und 4.4 
eine Liste der einzugebenden Konstanten). 


Bei den Programmbeschreibungen wurde davon ausgegangen, daß die Programme auf 
Magnetkarten gespeichert sind. 


Die Programmlisten haben folgende Form: 


000 43 RCL 
001 01 01 
002 42 STO 


003 02 02 


Dabei bedeuten die ersten drei Ziffern die Nummer der Programmspeicherstelle, die beiden 
Ziffern in der Mitte den Tastencode, rechts steht das Tastensymbol. 


0.1 Der Rechner Tl 59 


Der TI 59 ist ein programmierbarer Taschenrechner, der in drei Operationsarten betrieben 
werden kann: 


— Verwendung des Moduls 
— Durchführung von selbst eingegebenen Programmen 
— Benutzung als über die Tastatur zu bedienender Rechner 


Wichtig für den Leser ist die zweitgenannte Betriebsart. 
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Folgende Zeichnung zeigt die Speicherkapazität und -verteilung des Rechners 


88 
800 


| 


72 


480 


S 


Programmspeicher Datenspeicher 


Wird der Rechner eingeschaltet, stehen 60 Datenspeicher und 480 Programmspeicherstellen 
zur Verfügung. Mittels der Tastenfolge n 2nd Op 17 läßt sich die Speicherbereichsverteilung 
— etwa bei Programmen mit hohem Datenspeicherbedarf — auf 10°‘ n Datenspeicher ändern. 
Die zugehörige Anzahl von Programmspeicherstellen entnehme man jeweils der Skizze. Bei 
den einzelnen Programmen ist angegeben, ob und in welcher Weise die Speicherbereichsver- 


teilung geändert werden kann oder muß. 


0.2 Eingabe von Programmen 


Vor Eingabe eines Programmes empfiehlt es sich stets, den Rechner kurz auszuschalten, um 
sicherzustellen, daß alle Register gelöscht sind. Will man ein Programm in den TI 59 eingeben, 
geht man wie folgt vor: 


1. Einschalten des Rechners 

2. Taste LRN drücken. Es erscheint die Anzeige 000 00. 

3. Eingabe der in den Programmlisten abgedruckten Programmbefehle 
4. Taste LRN drücken 


Jetzt hat der Rechner das Programm gespeichert. 
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0.3 Magnetkarten 


In den Programmbeschreibungen wurde davon ausgegangen, daß die Programme von Magnet- 
karten in den Rechner eingelesen werden. Aus diesem Grunde und aus Bequemlichkeit ist es 
sinnvoll, ein in den Rechner eingetastetes Programm auf eine Magnetkarte zu übertragen. 
Eine Magnetkarte kann den Inhalt zweier Blöcke speichern; die Nummer der gespeicherten 
Blöcke sollte man beim Beschriften der Karte in den dafür vorgesehenen Kästchen links bzw. 
rechts oben auf der Karte vermerken. Das Beschreiben der Magnetkarte mit Block m durch 
den Rechner geschieht folgendermaßen: 


1. Tasten m 2nd Write drücken 
2. Magnetkarte einschieben 
3. In der Anzeige erscheint die Nummer des Blocks m 


Block m ist jetzt auf die Magnetkarte übertragen und wird mit jedem Einlesen im Rechner 
wieder abgespeichert. Das Einlesen einer Magnetkarte geschieht so: 


1. Anzeige auf O stellen 
2. Magnetkarte einschieben 
3. In der Anzeige erscheint die Nummer des eben eingelesenen Blocks 


Zu beachten ist, daß Magnetkarten nur dann eingelesen werden können, wenn die Speicher- 
bereichsverteilung dieselbe ist wie beim Beschreiben der Magnetkarte. 
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1 Matrizen 


1.1 Produktsumme 


Das Programm bestimmt die Produktsumme s = a'b zweier n-Spalten a und b. 


Programminstruktionen 


Tom TE [rm Tamm — 


Einlesen der Magnetkarte (Block 1) 


Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 4 


Eingabe von a,, ..., an; B1. ---, En 
Dabei bedeutet die Anzeige i, daß bisher 
i-1 Daten eingegeben wurden 


Ende der Koeffizienteneingabe 
Eingabe von k und n 


Ergebnisanzeige 


1.1 Produktsumme 5 


Beispiel 
3.1 4.2 
Man berechne s=a'b mit a=|2.7 | und b=|0.8|. 
1.5 2.3 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von: k 
a, 
a2 
a3 
b, 
b2 
b3 
Ende der Koeffizienteneingabe 


BNSOUPRODN A { 


Eingabe von: k 
n 


Anzeige von: s 


Programm 1.1 Produktsumme (Sl Produ I 


Ku 

ano TE LEL LE 215 a0 4 95 = 
act Al HA A ale 44 5 In Te BT 
lz ®1l Rei 699 [ir DE Eh [NE Sr u re 
a2 42 STD 97 DIE &1 5TOagk DE 01 A 
NE u Z Ip Sn ! z 4a 44 SLUM 
ads [ us so 00 co 
LE, 7 cSl 44 SUM 
alu} 1 LSEO DE GE 
Ale 4 055 47 IE 
cc 4 = = a4 01 01 
ao at 72 Rl# 055 13 
al a4) 00 G [SE 42 RÜL 
ciz 00 00, ndz 65 = [E77 DEE 
02 01 41 42 73 Rls 0a 91 Re 


144m) 44 DE 08 
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1.2 Matrizenprodukt 


Das Programm berechnet das Produkt C einer n,l-Matrix A = [a;;] mit einer I,m-Matrix 
B=[b;«]. © = [ci«] ist eine n,m-Matrix mit 


Cik = > ajj' bjk 


i=1 


Programminstruktionen 


Eingabe Anzeige 


Einlesen der Magnetkarte (Block 1} 
Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 10 


Eingabe der Koeffizienten zunächst von A 
spaltenweise, dann von B spaltenweise 


ni+1 
ni+2 


ni+Im+1 
Ende der Koeffizienteneingabe 0 


Eingabe von k, n, I,m k 
n 
| 


Ergebnisanzeige C1 


Registerinhalte 


Roo: ..-, Rog: Programmzeiger 

Rey... Re+ni-1: 9115 ---, Anl 
Resnts er, Res tnem)=1: 144 ++, Dim 
Re+itantm)ı Hr Rk+ı (ntm) +nm-1! ©1413 ++, Cnm 
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Beispiel 
5516|. 
Man berechne C=A-B mit A= 2 und B=|2 7 ı 
131 24 I 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von: K 

a1 

221 

412 


PIE) 
b23 
b33 
Ende der Koeffizienteneingabe 


0 
2 
1 
4 
3 
3 
1 
5 
2 
3 
1 
7 
4 
6 
1 
3 


Eingabe von: k 


_ Pe EEE EEE GE EEE NEE GEBE | 
ONDOO ON PD -OODODO NUT PWDND—-.-. 


N” 
S 


Anzeige von: 


DND#A- 
on» 


_ 
DS} 


: _I27 42 25 
Eric 11 26 a 
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Programm 1.2 Matrizenprodukt 
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2 Lineare Gleichungen und Ungleichungen 


2.1 Der Algorithmus von Gauß 


Das Programm berechnet die Lösung eines linearen Gleichungssystems Ax = a (A reguläre 
n,n-Matrix, a n-Spalte) nach dem Gaußalgorithmus. Dabei wird die Matrix [A, a] mittels 
geeigneter Zeilenumformungen in die Matrix [R, b] überführt. Die Lösung des linearen 
Gleichungssystems Rx = b, die durch Rückwärtseinsetzen gewonnen wird, löst auch das 
System Ax = a. Das Programm hält, wenn im Verlauf der Rechnung eines der Diagonal- 
elemente r;, von R zu Null wird, und zeigt dies durch eine blinkende Anzeige an. Vor dem 
Rückwärtseinsetzen gibt das Programm den Wert der Determinante von A aus. 


++ I 


Programminstruktionen 


2 


Magnetkarte einlesen (Biock 1, 2) 


Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 11 


Eingabe der Matrix [A, a] spaltenweise 


Ende der Koeffizienteneingabe 
Eingabe von k und n 


Anzeige von det A 
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(falls dt A# 0) R/S 
Ergebnisanzeige 


R/S 


Registerinhalte 


Roo. ---, Rıo: Programmzeiger 
Res... Rek+n2 1: 411, ---, Ann 
Re+n24 2. Rktnin+1)-1: a4, +, An 


Bemerkungen 


1. Das Programm „schreibt‘’ die gesuchte Matrix [R, b] über die eingegebene Matrix [A, a]. 

2. Bei gewöhnlicher Speicherbereichsverteilung bearbeitet das Programm Matrizen bis zur 
Ordnung n = 6; bei Änderung der Speicherbereichsverteilung auf 70 Datenspeicher mittels 
7 2nd Op 17 auch der Ordnung n = 7. 


Beispiel 

Gesucht ist die Lösung des linearen Gleichungssystems Ax = a mit 
2 20 6 

A=|2 1 Tt|und a=|7 
112 9 


Anmerkungen 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 
Eingabe von: k 


911 
271 
231 
412 
922 
932 
913 
923 
933 
a, 

a 


oSOD-O--.D-+-DDKND- 
ah ah ur 
OD-OOO Ss 1 PODODD—- DD 


a3 
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Ende der Koeffizienteneingabe 


Eingabe von: k 


Anzeige von: 


Anzeige von: 


Programm 2.1 Der Algorithmus von Gauß 


DL Bi =- NIE] NIR| 
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2.2 Der Gaußalgorithmus mit Pivotsuche 


Das Programm löst auch solche linearen Gleichungssysteme Ax = a (A n,n-Matrix, a n-Spalte), 
bei denen der gewöhnliche Gaußalgorithmus wegen einer Null in der Hauptdiagonalen der zu 
erzeugenden Matrix R versagen würde (siehe 2.1 „Der Algorithmus von Gauß’’). Beim Gauß- 
algorithmus mit Pivotsuche wird vor dem Eliminationsschritt j das betragsgrößte Element der 
„Restspalte’’ j 


a,;| = maxla; 
| rj! mex| jl 


gesucht und anschließend die Zeile r der Matrix [A, a] mit Zeile j vertauscht. 
Das Programm berechnet nicht den Wert der Determinante von A. Es bearbeitet Matrizen 
bis zur Ordnung n = 6. 


Programminstruktionen 


reisen _ [Em Tre] ame 


1 | Magnetkarte einlesen (Block 1,2) 
2| Programmbeginn 


3| Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 15 


4 | Eingabe der Matrix [A, a] spaltenweise 


5 | Ende der Koeffizienteneingabe 
6 | Eingabe von n und k 


7 | Ergebnisanzeige 


Registerinhalte 


Roo. -.-, Rıa: Prograrnmzeiger 
Res -.., Re+n2-1: 441, +, Ann 
Re+n2s ---, Rktn2+n—1: 91, :-:,@n 
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Beispiel 

Gesucht ist die Lösung des linearen Gleichungssystems Ax = a mit 
220 8 

A=|1 1 3| und a=|10|. 
132 14 


Magnetkarte einlesen (Block 1, 2) 


Programmbeginn 


_ 


Eingabe von: k 


411 
231 
831 
412 
422 
432 
413 
23 
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Die LR-Zerlegung ist ein konzentrierter Gaußalgorithmus, der die Matrix A zerlegt in das 
Produkt zweier Dreiecksmatrizen L und R. Dabei sei A zunächst ohne Zeilenvertauschungen 


zerlegbar, was i.a. nicht der Fall ist. 


Die LR-Zerlegung läßt sich zur Lösung linearer Gleichungssysteme verwenden. Besonders 
sinnvoll ist ihre Anwendung, wenn mehrere Systeme Ax =a,,..., Ax = a, mit identischer 
Koeffizientenmatrix A vorliegen. Zur Lösung dieser r Systeme wird die LR-Zerlegung ein- 


mal bereitgestellt. 


Das Programm wurde in zwei Teile zerlegt; dadurch läßt es sich auf Matrizen bis zur 


Ordnung n =6 anwenden. 
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Programminstruktionen 


Ce See 


„Teil 1° einlesen (Block 1) 
Programmbeginn „Teil 1" 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 11 


Eingabe der Matrix [A, a] spaltenweise 


Ende der Koeffizienteneingabe 
Eingabe von n und k 


„Teil 2” einlesen (Block 1, 2) 
Programmbeginn „Teil 2” 
Ergebnisanzeige 


"a 


Eingabe einer neuen „rechten Seite 


Ende der Eingabe 
Ergebnisanzeige 


Registerinhalte Bemerkungen 

Roo, ---, Rio: Programmzeiger 1. Ist eine LR-Zerlegung von A nicht möglich, so 
Res... Rk#n2-1: A445 ..-, Ann hält das Programm, und der Rechner zeigt dies 
Re+n2a.., Rk+n2+n—1: Ayser., An durch eine blinkende Anzeige an. 


2. Die Schritte 10 bis 12 können beliebig oft 
wiederholt werden. 


18 2 Lineare Gleichungen und Ungleichungen 


Beispiel 


Gesucht ist die Lösung der linearen Gleichungssysteme Ax=a und Ax=a’ mit 


2 2 0 6 12 
A=|2 1 1|, a=|7|, a’=|18 
1122 9 14 


„Teil 1’ einlesen (Block 1) 
Programmbeginn „Teil 1” 
Eingabe von: k 


911 
871 
231 
412 
422 
432 
913 
423 
433 
a] 

a, 

a3 

Ende der Koeffizienteneingabe 


OS OND-O 2..2.D- -ANDDND 


Eingabe von: n 
k 


„Teil 2’ einlesen (Block 1, 2) 
Programmbeginn „Teil 2’ 


Anzeige von: x4 
x2 
x3 

neue „rechte Seite” a’ 

Eingabe von: a} 
a2 

' 
a3 

Ende der Eingabe 

Anzeige von: x, 
x2 


x3 
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Programm 2.3 Die LR-Zerlegung 
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2.4 Die LR-Zerlegung mit Pivotsuche 


Die gewöhnliche LR-Zerlegung versagt, wenn eines der Diagonalelemente von R zu Null wird. 
Dann hat A keine LR-Zerlegung, wohl aber die Matrix P-A, wobei P die Permutationsmatrix 
der Zeilenvertauschungen ist. Das Programm bestimmt also Dreiecksmatrizen L und R mit 
L'R=P-A. Es wurde in drei Teile zerlegt und gestattet so die Anwendung auf Matrizen bis 
zur Ordnung n = 5. Für die Lösung mehrerer linearer Gleichungssysteme gilt das entsprechende 
wie beim Programm 2.3 „Die LR-Zerlegung”. 


Programminstruktionen 


7 Veran Ten [Tome] An] 


„Teil 1’ einlesen (Block 1) 
A 
k R/S 


Programmbeginn „Teil 1‘ 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 14 


4| Eingabe der Matrix [A, a] spaltenweise 


Ende der Koeffizienteneingabe 
Eingabe von n und k 


„Teil 2° einlesen (Block 1, 2) 
8| Programmbeginn „Teil 2” 


„Teil 3° einlesen (Block 1, 2) 
Programmbeginn „Teil 3’ 
Ergebnisanzeige 


Eingabe einer neuen „rechten Seite’ a 


Ende der Koeffizienteneingabe 
Ergebnisanzeige 


Registerinhalte 

Roo, -.., Rıa: Programmzeiger 
Re, ur Re+n2-1: 11, ---, Ann 
Re+n2, Pe Re+n24+n-1: 41, ..-, An 


Re+n2+#ns -, Rk+n?2+2n-ı: Zeilenindizes 


2 Lineare Gleichungen und Ungleichungen 


2.4 Die LR-Zerlegung mit Pivotsuche 


Bemerkungen 


1. In den Registern Rk+n2 +ns ---, Rk+n2+2n—ı werden die Zeilenvertauschungen 
gespeichert. 

2. Ist A nicht regulär, so hält das Programm und der Rechner zeigt dies durch eine 
blinkende Anzeige an. 


Beispiel 

Gesucht ist die Lösung der linearen Gleichungssysteme Ax=a und Ax=a’ mit 
220 6 12 

A=|1 1 2|, a=-/9|, a’=|14 
2171 7 18 


23 


„Teil 1” einlesen (Block 1) 


Programmbeginn „Teil 1" A 

Eingabe von: k 14 R/S 
a1 2 R/S 
231 1 R/S 
331 2 R/S 
a2 2 R/S 
422 1 R/S 
432 1 R/S 
413 0 R/S 
a23 2 R/S 
433 1 R/S 
a1 6 R/S 
Ep 9 R/S 
a3 7 R/S 

Ende der Koeffizienteneingabe B 

Eingabe von: n 3 R/S 
k 14 R/S 

„Teil 2° einlesen (Block 1, 2) 

Programmbeginn „Teil 2” c 

„Teil 3° einlesen (Block 1, 2) 

Programmbeginn „Teil 3" D 


Anzeige von: x; 
x2 R/S 
%3 R/S 


NO POND- -— 


© 


N —- 
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neue „rechte Seite’ a’ 


Eingabe von: aı 
Ep) 
a3 

Ende der Eingabe 

Anzeige von: x, 
x2 
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2.5 Inversion mit totaler Pivotsuche 


Das Programm bestimmt die Inverse A”' = [aj.] einer regulären n,n-Matrix A mittels des 
Austauschverfahrens. Dabei wird in jedem Schritt innerhalb der bisher nicht getauschten 
Zeilen und Spalten eine totale Pivotsuche durchgeführt. Der jeweilige Pivot wird in die Haupt- 
diagonale getauscht. Nach Durchführung der Inversion werden die Zeilen und Spalten in 
natürlicher Reihenfolge geordnet. 


Das Programm zerfällt in drei Teile und gestattet so die Anwendung auf Matrizen bis zur 
Ordnung n =6. 
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Programminstruktionen 


[Tem Teen Tran TA — 


„Teil 1” einlesen (Block 1) 
Programmbeginn „Teil 1” 
Eingabe der Matrixordnung n 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 11 


Eingabe der Matrix A spaltenweise 


Ende der Koeffizienteneingabe 


„Teil 2’ einlesen (Block 1, 2) 
Programmbeginn ‚Teil 2° 


„Teil 3° einlesen (Block 1, 2) 
Programmbeginn „Teil 3’ 


Ergebnisanzeige (A”' spaltenweise) 


Registerinhalte 


Roo, --., Rıo: Programmzeiger 
Re, ..., Re+2n-ı: Zeilen- und Spaltenindizes 


Re+2n: ang Rek+n2 +2n—1 © A411, +++, Ann 
Bemerkungen 
1. In den Registern Rx, ..., R«+2n-ı werden die Zeilen- und Spaltenvertauschungen 


gespeichert und für den Rücktausch von dort abgerufen. 
2. Ist A singulär, so hält das Programm während der Ausführung von „Teil 2’ und zeigt 
dieses durch eine blinkende Anzeige an. 


Beispiel 


Gesucht ist die Inverse der Matrix A = 


DB -anN 
a orN 
-nNo 
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„Teil 1 einlesen (Block 1) 


Programmbeginn „Teil 1" 


[ A 


Eingabe von: 
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Ende der Koeffizienteneingabe 


„Teil 2” einlesen (Block 1, 2) 
Programmbeginn „Teil 2” 


„Teil 3° einlesen (Block 1, 2) 
Programmbeginn „Teil 3” 


Anzeige von: a4, 


Esistalo A'=| 05 05 -1 
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Programm 2.5 Inversion mit totaler Pivotsuche 
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2.6 Die Cholesky-Zerlegung 


Eine reguläre, symmetrische, positiv-definite n,n-Matrix A läßt sich zerlegen in A=CTC, 
wobei C eine rechte obere Dreiecksmatrix ist. Wie die LR-Zerlegung verwendet man auch 

die Cholesky-Zerlegung zur Lösung linearer Gleichungssysteme Ax = a. Wegen der Symmetrie 
braucht man die Elemente von A unter der Hauptdiagonalen 
nicht abzuspeichern; eingegeben wird nur der in der Skizze 
schraffierte Teil der Matrix [A, a]. 


Das Programm wurde in zwei Teile zerlegt und gestattet so die 
Lösung linearer Gleichungssysteme bis zur Ordnung n = 8, bei 
Änderung der Speicherbereichsverteilung auf 70 Datenspeicher 
mittels der Tastenfolge 7 2nd Op 17 sogar der Ordnung n = 9. 


2.6 Die Cholesky-Zerlegung 


Programminstruktionen 


Eingabe 


„Teil 1'’ einlesen (Block 1, 2) 


Programmbeginn „Teil 1’ 


3) Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 11 

4| Eingabe der Elemente von A oberhalb und 
einschließlich der Hauptdiagonalen spalten- 
weise, anschließend Eingabe der Elemente 
von a 


5| Ende der Koeffizienteneingabe 


Eingabe von n und k 


„Teil 2° einlesen (Block 1, 2) 
8| Programmbeginn „Teil 2" 


Ergebnisanzeige 


Eingabe einer neuen „rechten Seite’’ a’ 


Ende der Koeffizienteneingabe 
Ergebnisanzeige 


Anzeige 


35 


=D 


IPB W@N 


n?+3n+2 


2 
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Registerinhalte 
Roo. ---, Rıo: Programmzeiger 
Rek,...,R n(n+1) : 811, 912, 922, 913, ---, dnn 
kt —— -1 
2 

R ntn+1ye rn, R n(n+3) "A157, @n 

gg En 
Bemerkung 


Ist eine Cholesky-Zerlegung von A nicht möglich, so hält das Programm und der Rechner 
zeigt dies durch eine blinkende Anzeige an. 


Beispiel 

Gesucht ist die Lösung der linearen Gleichungssysteme Ax=a und Ax =a’ mit 
210 5 1 

A=/11|4 1|, a=|9|, a’=|-1 
0 112 7 3 


„Teil 1” einlesen (Block 1, 2) 
Programmbeginn „Teil 1" 
Eingabe von: k 

471 

412 

432 

473 

23 

a33 

a, 

a 

a3 
Ende der Koeffizienteneingabe 


2 

2 
1 1 
2 2 
1 3 
4 4 
0 5 
1 6 
2 7 
5 8 
9 9 
7 0 

0 

0 


Eingabe von: n 
K 


nn © 


„Teil 2° einlesen (Block 1, 2) 
Programmbeginn „Teil 2’ 
Anzeige von: X; 

%2 

x3 
neue „rechte Seite’ a’ 
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u E 


a2 
az 
Programm 2.6 Die Cholesky-Zerlegung 


Ende der Koeffizienteneingabe 
Anzeige von: x 
x%2 
x3 
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Ü ı 1 07395 = 118.259 = 151 10 10 
03 EA 074 44 SUM 12 01 1 152 76 LEL 
I 2 STD 075 02 02 114 95 = 152 24 CE 
c or dr IrTE 97 188 115 42 5TO 154 43 REL 
038 25 CLR 077 09 09 11e 02 08 155 10 10 
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2.7 Die QR-Zerlegung und vermittelndes Ausgleichen 


Die OR-Zerlegung nach Householder führt eine n,m-Matrix Amit n>m und rang A=m 
durch Multiplikation mit einer orthonormalen Matrix Q in eine rechte obere Dreiecks- 
matrix R über. Der Algorithmus ist numerisch besonders günstig bei linearen Gleichungs- 
systemen Ax =a, deren Koeffizientenmatrix A schlecht konditioniert ist (d.h. deren Zeilen 
bzw. Spalten nahezu linear abhängig sind). Außerdem liefert die QR-Zerlegung die Lösung 
überbestimmter linearer Gleichungssysteme Ax = a, A n,m-Matrix, n > m, nach der 
Gaußschen Methode der kleinsten Quadrate (vermittelndes Ausgleichen). 


LH 


Das Programm wurde in zwei Teile zerlegt. In der folgenden Tabelle sind die bei normaler 
Speicherbereichsverteilung zulässigen Matrizenformate mit „+'' markiert; ändert man die 
Speicherbereichsverteilung mittels der Tastenfolge 8 2nd Op 17 auf 80 Datenspeicher, so 


sind auch die mit ‚‚o’’ gekennzeichneten Matrizenformate zulässig. 


N 


Programminstruktionen 


I Trrmen Ten Trans] Ama 


„Teil 1" einlesen (Block 1, 2) 2 
A 2 
k R/S 1 


2 |Programmbeginn „Teil 1" 


3 | Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 14 
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Teens _ [ersten] am 


4 | Eingabe der Matrix [A, a] spaltenweise 


mn+1 
mn+2 


Ende der Koeffizienteneingabe 
Eingabe von n,m, k 


„Teil 2° einlesen (Block 1) 
8 | Programmbeginn „Teil 2’ 


Ergebnisanzeige 


Registerinhalte 

Roo: .-., Rıa: Programmzeiger 

R«. OR Re+mn-1: A114, +++, Anm 

Rk+mn: si Re+mn+n—1 2 84, ++., dn 
Resmntns +, Rktmn+n+m-1: 81, +++, Sm 
Bemerkung 


Die $ı, ..., Sm sind Koeffizienten, die bei der Berechnung von R = QTA auftauchen. 


Beispiele 


1. Gesucht ist die Lösung des linearen Gleichungssystems Ax = a mit der nahezu singulären 
Matrix 
1113 14 
A=|I1 09 09| und a=|1 
0-01 0.2 0.3 
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„Teil 1’ einlesen (Block 1, 2) 
Programmbeginn „Teil 1” 
Eingabe von: K 

471 

221 

431 

412 

a22 

432 


ooONSOWP@QOD—- DD 


923 
433 
a] 
a2 


a3 
Ende der Koeffizienteneingabe 
Eingabe von: n 

m 

k 


„Teil 2° einlesen (Block 1) 
Programmbeginn „Teil 2” 


Anzeige von: xj 
%2 
x3 


2. Es soll eine Parabel y=co+cıt+t Cat? durch die vier Punkte der Tabelle 


t|-1 013 
vl 2123 


gelegt werden. Durch Einsetzen der Punkte in die Parabelgleichung erhält man das über- 
bestimmte lineare Gleichungssystem 


I ; 2 
100 u BE 
1119 a al Er 
124 e2 3 
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„Teil 1” einlesen (Block 1, 2) 
Programmbeginn „Teil 1’ 


_ 


Eingabe von: k 
471 
221 
a31 
841 
412 
422 
a32 
942 
413 
923 
933 
943 
a] 
a2 
a3 
a4 

Ende der Koeffizienteneingabe 


OD NP 2-09 2.2.2.2 0 DD 
ON PBODND-DDMN 


Eingabe von: 


„Teil 2’ einlesen (Block 1) 


Programmbeginn „Teil 2’ 


Anzeige von: Co 
cı 
c2 


Die gesuchte Ausgleichsparabel ist also 


y=13-0.1:t+0.5-t. 


44 


2 Lineare Gleichungen und Ungleichungen 


Programm 2.7 Die OR-Zerlegung 
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nr 


a 
COX 


43 
12 
05 
+ 43 RCL 
1 co du 
> 2 25 + 
= z ci 1 
: 2 5STD 2 s3 = 
ER] e LE 25 44 S11m 
234 r6 LEL 253 04 04 0505 
239 15 C' 254 44 SUM cl 1 
236 453 RÜL zs5 03 08 22 TNW 
gear ja 10 256 47 DSE a 
238 42 STO 257 09 09 
239 09 09 zes 17 EB!" 
240 Te LBEL 2592 43 RÜEL 
24 IT E" z£0 10 dl 


+ 
[&% 
ni 
= 


205 
Ace 


or 01 0 42 05 05 11# 
u0s ‚Ss > : 43 REL 113 
08 0 dl 01 01 114 
oo 54 5 73 - 115 
(dl ES x ci 1 116 
(lE 53 & 35 = 117 
113 43 RCL 42 53TD 118 
ol4 00 00 a 09 112 
DS 85 + TE LEL 120 
Sie 01 4 22 LN® 121 
Gr 54 43 RCL 122 
Dis 925 = 01 01 088 123 
01a 42 STO 75 - 033 01 01 124 
JEO 0303 42 RLL oa 75 - 125 
21 85 + 0309 131 43 REL 126 
(22 43 RÜL ss = 032 09 09 127 
23 00 00 42 STD 0192 95 = 128 
95 = us 08 024 44 SUM 123 

42 5T0 FE LEL 035 040g 130 

04 04 z2 IMY DEE 65 x 131 

73 RC 73 Rl# D" 43 PEL 132 

03 03 02 03 033 00 00 133 

22 IHM“ 85 x 033 75 - 134 

54 FI* 73 Ro 100 01 1 135 

04 04 4 04 101 95 = 1:36 

E 137 
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2.8 Zyklische Relaxation 


Zu einem linearen Gleichungssystem Ax =a (A n,n-Matrix) sei eine Näherungslösung p mit 
dem Residuum 


r:=Ap-a#0 


gegeben. Die Idee der Koordinatenrelaxation besteht darin, eine Komponente P; der Nähe- 
rungslösung p so zu ändern, daß die zugehörige Komponente r; des Residuums r zu Null 
wird. Bei diesem Algorithmus werden alle Koordinaten p, der Reihe nach so oft abgearbeitet, 
bis die Tschebyscheff-Norm des Residuums |Ir||_, eine vorgegebene Toleranz 0 >O unter- 
schreitet. 

Das Programm bearbeitet lineare Gleichungssysteme bis zur Ordnung n = 5, bei Änderung 
der Speicherbereichsverteilung auf 80 Datenspeicher mittels der Tastenfolge 8 2nd Op 17 

bis zur Ordnung n = 7. Si 


Programminstruktionen 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 10 


Eingabe der Matrix [A, a, p] spaltenweise 


n?+1 
n?+2 


n?+n+1 
n?+n+2 


n?+2n+1 
5 | Ende der Koeffizienteneingabe 
6 | Eingabe von o,n,k 


7 | Ergebnisanzeige 
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Registerinhalte 


Roo, ---, Roo: Programmzeiger Re+n2+ns +, Rk+n2+2n-1: P4r + Pn 

Re, u Re+n2-1: 411, ..-, Ann Rek+n2+2n, 4 Rk+n2+3n-1: My. in 

Re+n2 ER Re+n2+n—1: 941, .:+, An 

Beispiel 5 :4>0 2 

Zu dem linearen Gleichungssystem Ax=a mit A=|1 4 1| und a=|B8| seidie 
01 2 

Näherungslösung p=| 1.2 | gegeben. Gesucht ist eine verbesserte Näherung x mit 


-0.7 
Ir, = IAX-all_ < 0.01. 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 


_ 


ODD-O 2-0 -N0O 


Eingabe von: K 
271 
271 
231 
412 
@22 
a32 
413 
433 
433 
a, 
a 


ON OO UI PR@OND— DM 


a3 
Pı 
Pa 
Pa 
Ende der Koeffizienteneingabe 


Eingabe von: 0 


-0.00234375 
2.001171875 
-.0005859375 


Anzeige von: 


xı xıxı = 3 
ND 


0 
Die exakte Lösung ist x =| 2 
0 
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Programm 2.8 Zyklische Relaxation 
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DIR E] T = £ 
ac be [7 7 


a1l 
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cc cn 


n r3 Rs 
ul: u u u Eu Re 
cd 3 25 = 

5 r4 5M# 
Re u 


.cı 


11 


oe 


Mr) £ı 
TEN La Pr 


1 
1 
1 
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1 
123 
1 
1 
1 
1 
1 
1 
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- 
BL 
Di Br DE 
{= 
o 
Zı 


- 

un 

Zı 
(VUN: 
(BR ] 

m 

Der] 

m 


151 oo oc 44 Su 

152 925 = 105 os 
152 42 5STD EER + 

124 04 04 3 42 FÜEL 
153985 + + 00 00 
156 43 ROL ee - 

157 00 oc E44 SUM 
155 935 = r 0208 
152 42 5TO SE 752 
le 05 05 03 09 


DEU) 


fr 1 


or 
<ı 
.L 


03 

ler? TE LEL 

250 185 19 D" 
us 1E9 73 Rls 
IS2 43 111 & 17005 05 
053 00 00 112 42 171 34 #7- 
1S4 95 = 112 02 02 lr2 59 * 
135 42 5STD 114 85 + 173 73 Rl* 
SE 04 (d 112 02 5 ir4 02 02 
15776 LEL lilE 85 x 175 935 = 
13 C 117 43 RÜL ie 42 ST0 
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44 Su 342 RÜL 4 ug O4 04 
05 on ac +2 FÜL 1 Ess 
ar EHE INT In ci 1 
Duke] SS + 42 STD 44 SUM 
18 43 FÜL 03 09 a4 04 
42 Ts In) F& LEL Sr 15E 
ci gg = er 093 09 
5 42 5STD r2 FRü# ir BE" 
1 Rss 


2.9 Methode des verstärkten Abstiegs 


Diese Methode zur Lösung eines linearen Gleichungssystems Ax = a mit symmetrischer und 
positiv-definiter n,n-Matrix A ist ein Relaxationsverfahren (siehe 2.8 ‚„Zyklische Relaxation’), 
bei dem die Näherungslösung p nicht koordinatenweise, sondern in Richtung des Residuen- 
vektors r= Ap-a geändert wird. Der Algorithmus endet, wenn die euklidische Norm des 
Residuums r eine vorgegebene Toleranz 0 >O unterschreitet. 

Das Programm gestattet die Bearbeitung von linearen Gleichungssystemen bis zur Ordnung 
n=5, bei Änderung der Speicherbereichsverteilung auf 70 Datenspeicher mittels der Tasten- 
folge 7 2nd Op 17 auch der Ordnung n =6. 


Programminstruktionen 


[Team _ Tem ran] Aue 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 10 


Eingabe der Matrix [A, a, p] spaltenweise 


n?-+n+1 
n?+n+2 


n?+2n+1 
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Tremmn _ [em [ran T Amis — 


5 | Ende der Koeffizienteneingabe 


6 | Eingabe von o,n,k 


7 \ Ergebnisanzeige 


Registerinhalte 

Roo, -.-, Rog: Programmzeiger 

Res +-., Rk+n2-1: 441, .-+, dan 

Rek+n2 4... Rk+n2+n—1: 91, .-«, An 
Rek+n2#nı +, Rk+n2+2n-1: Pır::Pn 
Re+n2+2n, ++, Rk+n2+3n-1: Ms +--, In 
Beispiel 


210 2 
Zu dem linearen Gleichungssystem Ax=a mit A=|1 4 1| und a=|8| sei die 
0 12 


0.5 
Näherungslösung p =! 1.2 | gegeben. Gesucht ist eine verbesserte Näherung x mit 
-0.7 


IIrll, = IIAX - all, < 0.005. 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 


_ 


D-O-.2-.0->D1D0O 


Eingabe von: k 
a1 
821 


431 
412 
922 
932 
913 
423 
433 


SOC GT NO TI PRPOND-DNMN 


— 
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a4 2 


11 
42 12 
a3 


P; 
Pa 
Pa 


Ende der Koeffizienteneingabe 
Eingabe von: o 


Anzeige von: .0009332064 
1.998317586 
.0009243247 


0 
Die exakte Lösung ist x =| 2 
0 


Programm 2.9 Methode des stärksten Abstiegs 


[NIRTE| ? 


r 


I" tm 
r 


a 


DAY er u Dep) 


Ce BE ui 
DO in 


003 
old 
BB 
N12 
12 
214 
015 
Üle 
IT 


u eu 
Mn 
(Be I 


m+ m 


1 5: 
1 RS 0S4 0 ci 44 +5 
42 STD Dun] ss 4% sl 72 ST8 
d2 02 STD Se 0 us2 07 07 
ci 1 IE [Ki r 0Ss3 01 1 
42 STD CLR Iubrt=} DE4 Sl 
0302 R=3 usa 7 nen 
4: STD Te Br) DEE 
[1 & 
D & 
[et 


a 


3 
& 
S 
e) 
5 
=) 


El 


rs] 


17 
De We 
mr 
nn ı 5 SE 
Tl 
ra 
u! 


og4c Du 
ü4l 32 
4 35 
a4: 425 


DiY} 
un u | u LE) 
m 


] 
DEN u FE = 
vn 
+ 


Ir 

5 
ı 
r 


Ib Er ü44 07 ou 44 SUM 
211900 ü4S 553 + 6 6 
NS Er = Idee 43 RÜL O1 1 
1 76 47° 000 0 44 ZUM 
13 üdE 95 = TS TE 
25 ü43 42 3T0 Ar DSE 
Ep usa de 08 (2003 
SE Sl Te LEL a4 75 
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104 00 cu 42 FÜL 


110 08 06 

111 85 + 

112 43 FÜL 

1132 00 00 

1id 23 8 

113 95 = 

ile 42 5T0 

ir or 07 

115 85 + 

11a 0202 

120 685 8 

121 43 REL 

122 00 Gü 

123 95 = 

124 42 570 

125 08 08 Er 

126 75 LEL a = 

le? 35 178 42 5STD 

1ES 43 REL Tu NT} 

129° 00 0 35 + 

130 42 5TO 42 REEL 

1231 03 03 Do 00 

132 00 0 us = 

133 r2ST# 42 ST0 

134 08 08 as 02 

1355 76 LEL FE LEL 

126 42 REL 44 SLUM 

138 1]-} LE 7 7 + 

13? 683 8 ES x url 

140 73 Rl* 42 FOL 32 mE 
141 07 07 03 09 SI + 
142 95 = us = 43 FÜL 
1432 74 SM& 74 SH 1 01 
144 08 08 lE 0% gs + 
145 43 RCL 72 Re 42 FELL 
ide 00 00 02 02 00 00 
14” 44 Sum Ba“ u = 
1482 068 06 42 RÜL 42 5STD 
143 01 1 0309 Ü [3 
150 44 5UM 35 = re LEL 
151 07 07 74 SMs Ss2 EE 
152 97 DSZ v2 Rlc# 


DE a N u 
ala a ge Pa ed Pl 
- 
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2.10 Lineare Optimierung 


Das Programm berechnet die Lösung eines linearen Programms in Normalform 


y} >20 
y„"By,+b >20 
z=c'y\,tc-max 


nach dem Simplexverfahren. Dabei sind y, und c m-Spalten, y, und b n-Spalten und B ist 
eine n,m-Matrix. Das lineare Programm wird spaltenweise als eine Matrix L der folgenden 


Form eingegeben: 
B | 
Derflı 


1- 


nn 


Das Programm wurde in drei Teile zerlegt. Es gestattet die Lösung solcher Optimierungs- 
probleme, deren Formate in der folgenden Tabelle mit ‚‚+‘’ gekennzeichnet sind. 


2 
3 
4 
5 
6 
7 
8 
9 


+ ++++++ + 
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Programminstruktionen 


„Teil 1” einlesen (Block 1) 
Programmbeginn „Teil 1’ 
Eingabe von n und m 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 13 


B|b 
Eingabe der Matrix L = Ed spaltenweise 
e'|c 


m{n+) 
m{n+1)+1 
m({n+1)+2 


{m+1) (n+1) 


„Teil 2° einlesen (Block 1, 2) 
Programmbeginn „Teil 2" 


„Teil 3° einlesen (Block 1) 
Ausgabe der primalen Lösung y, 
Ergebnisanzeige 


Registerinhalte 


Roo, --., Rı2: Programmzeiger 
Rx, ..., Re+m+n-ı: Zeilen- und Spaltenindizes 


B|b 
Ret+m+n: »--, Rk+2m+2n+mn: Koeffizienten von L = ER 
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Bemerkungen 


1. Ist man nicht an der Lösung des primalen, sondern des dualen linearen Programms 


V2 <o0 
vi =-B'wt+te <o 
w =-b'wy+tc>min 


interessiert, so wird Schritt 9 der Programminstruktionen ersetzt durch Schritt 9a 


| 9a | Ausgebeceraustentum | — — ÄÜPD| | 


B|b 
Auch in diesem Fall wird die Matrix L = | = eingegeben. 
c c 


2. Existiert keine optimale Lösung, weil die Menge der zulässigen Lösungen unbeschränkt ist, 
hält das Programm und der Rechner zeigt dies durch eine blinkende Anzeige an. 


Beispiel 


Gegeben ist das lineare Optimierungsproblem 


10 
Yı|l + 124 20, z=2y,ty,r0>ma; 
Y2 6 


„Teil 1’ einlesen (Block 1) 
Programmbeginn „Teil 1" 
Eingabe von: 
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2 
3 
4 
5 
6 
7 
8 
9 


„Teil 2’ einlesen (Block 1, 2) 
Programmbeginn „Teil 2” 


„Teil 3’ einlesen (Block 1} 
Ausgabe der primalen Lösung 
Anzeige von: Zmax 

Yı 

Y2 


Programm 210 


DC 76 LBL 44 01 1 DEE 085 
cool dl A 045 44 SuM ler nl 
Doz 21 Res Ode 05 05 ES 54 
03 42 5TO 047 44 SUM IT 
ca4 00 cc od 08 08 Dre Br Sc 
005  4= 5TO 43 97 DSE 071 04 
DT u Pe u ne 050 02 WE ürz2 TE 
u07 91 R=S 22 LN® Or 24 
Ice 42 5TO O1 1 Pi 42 
03 01 0 42 STD ? 

oil 42 STD 02 03 

Oli 04 04 on u 

012 921 RzS 53% 

O12 42 570 43 FÜL 

ül4 0202 on au 

015 42 5T0O E59 + 

Die 05 05 O1 1 

u Er 1 u 54 5 

O1lS 42 5STD 63958 

012 06 DE 52 « 

Den TE LEL 43 FEL 


021 22 IM 
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"E LEL 0SE 85 x 1l&E 87 IFF 174 05 05 
u 052 42 REL 117 08 06 175935 = 
42 ROL 0&0 01 11& 42 5TD 176 42 5T0 
oil 01 os1 112 42 STO Ir7 03 03 
+2 5STD 062 120 08 08 lt 75 - 
a4 04 082 121 Z2 IM 172 43 RÜEL 
SS + 064 122 8£& STF 120 05 05 
[u oe 065 22 06 08 181 55 + 
ES ÜBE 124 43 RCL 182 43 RÜL 
42 RÜL 135 02 09 123 01 01 
u 5) 126 42 STD 184 85 + 
I + 127 5 185 42 FEL 
42 RCL 128 61 GTD 1E& 06 
02 08 [ 123 235 1°X 127 95 = 
3 = 3x 120 76 LEL 122 42 570 
42 5STD 42 RÜL 121 42 ST0 123 04 04 
2 02 ar4 0505 43 RCL 120 73 Ro 
23 cp U7S 85 + 132 10 10 1231 02 02 
I 076 42 RÜL 134 75 - 132 63 EX# 
42 5TO 07 05 05 125 42 RCL 123 04 04 
us 0S Ir 73 - 126 08 085 124 63 Ex# 
Fe LEL ur od 1 13795 = 135 03 03 
25 CLR [u 15 5 = 1328 7 GE 136 42 RÜL 
r3 Rl* 081 42 5T0 139 35 1°X 137 00 m 
j 03 03 0S2 04 04 140 43 RCL 135 85 + 
25 94 +#- oe3 01 1 idi 10 10 122 01 1 
ZEO BEIM us 42 STD 142 42 STO zu0 95 = 
DE? FT? GE 0909 id3 082 08 2 42 STD 
Ü2S 32 Nr 43 RÜL 144 43 RCE 2 3 02 
023 01 1 ou 00 145 093 09 z 43 RÜL 
020 44 SUM 42 5STD 146 42 3T0 2 O2 08 
31 05 05 ur 07 147° 06 06 2 33 + 
32 85 + 00 0 148 76 LEL 2 43 FÜL 
33 43 RÜL 42 149 35 17% r Oi lt 
0 00 a5 1SDd 0 1 2 5 + 
935 = 26 151 44 SUM 2 42 RÜL 
44 ZUM 0% 152 03 03 os 05 
03 032 F& 153 44 SUM EI H 
38 97 137 33 154 04 04 43 RÜL 
039 04 04 73 155 44 SUM 00 0 
O4ü 25 CLR ud 155 09% 09 85 + 
O4l 25 CLR 094 22 IMW 15797 152 43 FELL 
042 91 RS 100 77 6 153 07 07 05 05 
042 TE LEL 101 34 TH 159 23 x8 #5 - 
44 32 X 102 81 GTO 160 232 IMY od 
45 43 RCL 10325 1% iEl 87 IFF 25 = 
Ode 02 08 104 FE LEL 182 0& 06 42 STD 
047 85 + 105 24 7 1E3 43 RCL 0309 
c4E 43 RCL 10E FE RC& 1E4 81 GTO 55 + 
043 00 ac 107 03 03 1865 16 RA! 43 FCL 
0so 85 + 108 55 = 1ER 7E LEL 06 0% 
051 43 RÜL 103 73 Rc# 167 42 RCL 75 - 
0S2 00 00 110 04 04 1E2 42 RÜCL ol 1 
053 65 8 1il 95 = 3 02 0 25 = 
054 42 RCL 11E 94 +7- 1r0 75 - 42 STD 
055 01 01 113 2 STD 171 01 1 U7 Fr 
056 85 + 114 10 10 lrE 85 + r3 Rl# 
0S7 02 2 115 22 IM 173 43 RCL Dr 07 
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KEN a Tr PRINT 
m 
Ina} 
ı ine m) m ( 


42 5STD TE LEL 42 FÜL 4153 75 - 

O7 07 S2 EE DE 08 u 1 5 u | 

O1 1 42 FEL 3 > 25 = 

r5 - 12 12 cl 1 42 STD 
42 RÜL Er o Eu „3 = 0309 
Ü DE Da +2 5T0O ci 1 

3 = 2 Rl& 10 10 Fe 

42 5STD 0308 I 42 FÜL 
1 11 TE rd - MT NT 
Y 'L > Rl# 42 RÜL sr 

44 SLM a 10 (5 05 = 5STD 
4 + a = 0 10 
I 11 REL 42 5T0 LEL 
£ Eu in IT 0309 5TO 
4 5 = re LEL RÜL 
43 RÜL - 45 


Ju 


[1 
Zen 
Arte gmoje gig 


Ta en Lin 
m 
zZ 
DER ER ı a Br BE a u 2 
rm 
Dee} 


SE Es ea 


55 + [1 Ü 13 F 

Ib Es LEL E Ei RÜL 
ge = [1% S ri urn AT 
42 5STD 43 ECL 423 RCL zZ IMW 


11T ll 


CE DE ua EN RRN 
+ 
ı 


= 
= 
o 
fer 
< 


$ 

+ 

LN [47] 
GC ı 

PR, 


DZ 


TEN Le 
+ “Lı 
r in 
Dir} 
Z 


xı 
N ur BR ı Ba Wera. Du OUETTE 


zZ 
je 
oe m 
IL Fre IN IT 
{1} 


2 + 44 SLM a3 09 Ar DSE 
42 FELL 12: 12 Be} + 03 03 
0& 06 37 D5Z 43 FEL &1 GTO 
Jg = 4 04 DT u ln) 423 ELCL 
4z 5T0O SZ EE 35 = ÜE DZ 
ic ic re LEL 29 + 
43 FÜL 43 78 43 REL 
O2 DE oil 1 ci 01 
E59 + 44 SM EI + 
43 RÜL 11 11 42 RÜL 
1 cl 44 Sum 00 00 
S + 3 03 5 Kr 


Zur 


2 2 97 082 +2 RCL 
EEK 03 03 2505 
43 RCL 44 SUM 33 + 
on 42 FÜL >33 +3 RCL 
= - oc m +4oc 0505 
43 RCL 25 + 401 ae 
03 03 1 1 402 2 2 
En + 35 = 42 E59 + 
Oi 1 2 5STD 404 42 REL 
35 = 203 405 DE NE 

2 ST 42 RCL 406 85 + 33 = 

02 02 407 42 RCL 42 STO 
ES + goes 00 Don 08 08 
43 RÜL 409 65 x 43 RLL 
Oi 01 410 43 FÜL ur 07 
EI + 411 05 05 35 1. 
42 FÜL dizE 89 + Y2ST# 
NT Te 413 43 FÜL 5 08 
35 + 414 05 05 12 EB 
29 CLR 
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ad FE LEL 0SzE 42 FÜL 104 65 126 TS D+ 
cool 12 c 053 04 04 103 92 4 iS7 00 0 
OO 42 FRCL 54 25 = Ile 43 FÜL 135 91 R#S 
03 02 08 üss £&7 E08 ir 01 152 &1 GTO 
Da4 85 + SE 21 RST 10385 + iel0 Fü RAD 
DoS5 085 os? ot 1 102 42 REL 181 TE LEL 
DOES NS5 44 ZUM 110 00 00 IEE TOR 
007 53 & us? 05 05 111 54 5 163 2 REEL 
DE 43 RÜL Del 37 1052 112 8275 + led 02 DE 
003 01 01 it 0606 113 42 RCL 165 35 + 
Dlü 85 + [ 73 - 114 00 00 166 53 %£ 
coll 43 RCL on u 115 65 1er 01 1 
colzE 0c oc I lle 423 FRCL lES 85 + 
13 54 0 &1 GTO 117 01 0 169 43 RÜL 
25 + ES + 1liö 95 = 170 05 05 
42 FOL TE LEL 11? 42 STD 171 75 - 
Dee IR) 21 RST 120 0707 172 42 FEL 
BE ® 43 RÜL 121.» 23 173 02 082 
43 RLL us 05 lzEZ ir4 54 
or i 25 + 122 175 865 x 
ag = 2 43 RÜL 124 irE 53 ı& 
Pr orz2 00 Du x 1FfF 43 RÜL 
ur 074 585 + 2 175 00 00 
"2 07543 REEL g ira 85 + 
7 TE 1 zZ 150 01 1 
Er 7 = = is1i 54 ı 
43 2 182 85 + 
[1 ] \ 1823 42 REL 
42 oc 1 124 00 00 
Inbe 031 133 185 85 + 
nBl 082 124 186 42 RÜL 
42 23 135 157 01 01 
ü4 0s4 136 is8 75 - 
FE ugs 137 152 01 1 
rt OSE 135 120 95 = 
+3 DEr 123 = 131 42 STD 
RS as 140 , - 
35 ht) 141 LEL 
42 Dan 142 z+ 
033 01 031 143 Rö# 
ode 95 = 092 144 


c4l 42 5T0 132 
de 05 05 cf 
oda 42 FÜL 1135 
a44 00 0 USE 
ü4S 0 42 5T0 37 
Ode 06 0E 198 
"47 Te LEL 038 
24 
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050 


Si 
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TE 
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44 
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?1 SER 148 
Er! 
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42 R 
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35 


Ba Be 5 
[a 2 | 
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pe 
Dee] 
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ı 
Fu 
Zı 
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Den 
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3 Iteration 


3.1 Vektoriteration nach von Mises 


Verfügt man über einen geeigneten Startvektor y, und besitzt die n,n-Matrix A einen betrags- 
größten Eigenwert X,, so konvergiert die Iterationsfolge 


1 
VeRYHFTTTH 12 
j le 
gegen den Eigenvektor x; von A und die Folge 
T 
Y;Y; i 
T > Wels 
Y; Yı-ı 


gegen den Eigenwert X}. Das Programm bricht ab, wenn 
Iy; y_ı' Iv;I_l. <e 


ist (€ >O Toleranz) oder die vorzugebende Maximalzahl N von Iterationen durchgeführt 
worden ist. Es bearbeitet Matrizen bis zur Ordnung n = 6. 


Programminstruktionen 


Term Tee Tre m — 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 11 


Eingabe der Matrix [A, y,] spaltenweise 
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Tees Terme [rm] ame — 


5| Ende der Koeffizienteneingabe 


6 | Eingabe von n,k, e, N 


7| Anzeige von X} und x, 


Registerinhalte 

Roo, ---, Rıo: Programmzeiger 

Res -.., Rk+n2-1: 911, ++, Ann 

Re+n2 ++, Rk+n2+n-1: Yors a 
Resnz#nı + Rk+n2+2n—1: yin LE 
Beispiel 


Mit höchstens N = 5 Iterationsschritten, der Toleranz e = 0.1 und dem Startvektor 
Yo=[1,0, 0]" soll der betragsgrößte Eigenwert X, und der zugehörige Eigenvektor x; 
der Matrix 


32 -1 
A=|26 -2 
00 2 


näherungsweise bestimmt werden. 


Magnetkarte einlesen (Block 1, 2) 


Programmbeginn 

Eingabe von: k 
911 
871 
331 
412 
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32 


Ende der Koeffizienteneingabe 


Eingabe von: n 


Anzeige von: X; 6.999746256 


5047690015 
1 
0 


Die exakte Lösung ist A}= 7 und x; = | 1 


Programm 3.1 Vektoriteration nach von Mises 


000 Te LEL o13 00 00 42 STO US? 73 Ro# 
oa 11 A 0208 02 06 uss 02 02 
002 91 Res 7& LBL 040 25 CLR 059 42 510 
002 42 5T0 iz 6 ö4l 91 ReS sc 10 10 
004 00 00 25 CLR o42 42 STO Osl 76 LEL 
os co 1 a1 043 05 05 DEE 32 INW 
ocs 42 STO 42 35T O44 42 RÜL sc | 

or 0 00 00 45 00 0 E44 BLUM 
002 43 RCL 5 - de 32 ME "65 02 02 
003 01 01 1 04T 85 + E86 73 Rl* 
010 1 Res = 042 42 RÜL 067 02 02 
Dll FE ST# I STD ö43 01 01 068 50 II 
SIE 00 0 1 03 050 95 = 069 22 IMW 
a2 01 1 032 LR: 051 42 STD Or0 7T GE 
Old 44 SUM 033 15 us2 02 02 Ürl 23 INH 
015 00 00 034 053 72 Rl# ürz 67 EU 
Ole 44 SuM 035 054 02 02 073 23 LNX 
17 01 m 26 055 50 I=]l 074 32 MIT 
0lS 61 GTO NET USE 32 HIT 75 73 Rl#z 
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& 42 570 5 12 73 $ 
rF 02 08 Be 13 02 > 
g 59 130 0 14 33% 
a De 15 44 Sum 
0 00 192 85 + 1 08 ds 
1 334 123 42 RÜL - EEE T 17 73 RC 
e 95 = 194 01 ; 6 LE 18 02 02 
3 42 STO 135 = >37 17 B a - 
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5 99 + 137 In Pal >53 ?=3 Rl# n 03 03 
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a7 DS 42 FCL 

09 09 oa 00 

16 RA’ 42 53T0 

325 42 RÜL 0309 

339 02 08 ZZ IHW 

a40 55 = 44 51M 

41 43 FÜL E02 
242 ur 07 42 RÜL 3 
242 95 = 10 10 03 
244 42 570 1 Re ER 
245 10 10 TE LEL e 


3.2 Inverse Iteration 


Das Programm berechnet den betragskleinsten Eigenwert X, der n,n-Matrix A als betrags- 
größten Eigenwert N von A"', Dabei wird die Iterationsvorschrift 
= A'y_, 

ersetzt durch das Lösen des linearen Gleichungssystems 

LRy;=y_,. 
wobei die LR-Zerlegung von A vorher durch das Programm 2.4 „Die LR-Zerlegung mit 
Pivotsuche” bereitgestellt wird. Das Programm bricht die Iteration ab, wenn der Wert 

Iv; TY=4 lv; 


eine vorzugebende Toleranz e >O unterschreitet oder die Höchstzahl N von Schritten durch- 
geführt wurde. Es ist in drei Teile zerlegt und auf Matrizen bis zur Ordnung n = 5 anwendbar. 
Damit die Speicherbelegung übereinstimmt, muß beim Programm „Die LR-Zerlegung mit 
Pivotsuche’’ die Nummer des ersten zu belegenden Speicherplatzes k = 14 gewählt werden. 


Programminstruktionen 


Tee Tem ra] ame — 


„LR-Zerlegung mit Pivotsuche — Teil 1” 
einlesen (Block 1) 


Programmbeginn „Teil 1'' 
Eingabe der Nummer des ersten zu 


belegenden Speicherplatzes k = 14 


Eingabe von A spaltenweise 
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5| Ende der Koeffizienteneingabe 
6| Eingabe von n und k 


„LR-Zerlegung mit Pivotsuche — Teil 2” 
einlesen (Block 1, 2) 


Programmbeginn „Teil 2’ 


„Inverse Iteration — Teil 1’' einlesen 
(Block 1) 


Programmbeginn „Teil 1” 


Eingabe des Startvektors Yo 


Ende der Koeffizienteneingabe 
Eingabe von N und e 


„Inverse Iteration — Teil 2° einlesen 
(Block 1, 2) 


Programmbeginn „Teil 2 


„Inverse Iteration — Teil 3° einlesen 
(Block 1) 


Programmbeginn „Teil 3’ 


Anzeige des Eigenwerts X, und des 
Eigenvektors x 


Registerinhalte 


Roo, ..., Rı3: Programmzeiger 
Rı4, ---., Rn2+13: Koeffizienten ee LR-Zerlegung von A 
Rn2+14,.-., Rn2+n+13: yü ven yÜ 


Rn2 +n+14: +++, Rn2+2n+13: Zeilenindizes 


1 1 
Rn2+2n+14, ---, Rn2+3n+13: yir nn ee ! 


Rn2 +3n+14, ---, Rn2+an+ 13: Zeilenindizes 


66 3 Iteration 


Beispiel gt 

Gesucht ist der betragskleinste Eigenwert Aa der Matrix A=|0O 4 2 | und der zuge- 
123 
e=0. 


hörige Eigenvektor x. Dabei sei yg=[1,0, 0], N =10 und 


„LR-Zerlegung mit Pivotsuche — Teil 1” 
einlesen (Block 1) 
Programmbeginn „Teil 1’ 
Eingabe von: k 

947 

21 

431 

a2 

422 

432 

413 

423 

433 
Ende der Koeffizienteneingabe 


OD-DPO-0O0—-.>» 
OOO NO TNIT PO DD 


_ 


Eingabe von: 


„LR-Zerlegung mit Pivotsuche — Teil 2'' 
einlesen (Block 1, 2) 


Programmbeginn „Teil 2” 


„Inverse Iteration — Teil 1’ einlesen (Block 1) 


Programmbeginn „Teil 1’ 


Eingabe von: yP 


yo 
yo 


Ende der Koeffizienteneingabe 


OO 2099 N 2-60 


_ 


Eingabe von: N 
€ 


oO 


3.2 Inverse Iteration 67 


„Inverse Iteration — Teil 2” einlesen 
(Block 1, 2) 


Programmbeginn „Teil 2” 


„‚Inverse Iteration — Teil 3’ einlesen (Block 1) 

Programmbeginn „Teil 3'' 

Anzeige von: Aa ‚354248689 
xı 1 
X2 ‚3542448811 
Xa -.6457466823 


Programm 32 


NINE 
IIuBI 
Oz 
002 
nun‘ 
[nIRt=] 
le 
our 
ce 
207 
alu 
call 
clE 
013 
als 
(15 ! | 
üle FE LEL 
olr © 
J1lE 
19 
np e) 


n 


+ 


n 
Inne 


Eu 


ızı 
(Bee We Bun 
[3 EN ı E DEE 2 


De NE 


u 


ıı 
(ei 


nie "oo Del 


me 
Dee) 


Dol 786 LEL a0s ie A" üle 95 = 
ool 13 C [alu )=) 2 5STD d1lr 42 
002 28 CF 0i0 05 05 ls 07 
003 43 REL oil 85 + 012 76 
oco4 00 00 alE 02 2 0 233 
005 42 5STO 012 65 8% r3 
le 04 04 o14 3 FÜL 05 
o07 71 SER 015 00 53 + 


68 


032 
033 
034 
035 
USE 
037 
28 
{129 
40 
041 
ü42 
043 
044 
045 
046 
047 
048 
049 
050 
051 
052 
053 
054 
055 
056 
057 
058 
059 
Deo 
061 
062 
063 
064 
065 
066 
067 
068 
069 
070 
or 
072 
073 
074 
075 
076 
077 
078 
079 
080 
081 
082 
083 
084 
085 
086 
087 
088 
0893 
020 
091 
092 
093 


SUM 
05 
SLIM 


034 
0139 
(36 
üsr 
098 
094 
Loc 
101 
108 
103 


a ı Bo SSL DO 


Pal 2 Cd CD LA LOG Da PA Tu Pi TI Pe PP Pc Pe me em mn nn m 
ra IN Ba Dre TI TIN Be c 


ba a a na 


2u3 
204 
205 
206 
207 
208 
209 
210 
zii 
212 
213 
214 
215 
216 


21? 


roTaro to oa 
MAN sl Zueuc cc 


fa 3 Ri To Fo fa Pa Po Pf 


LI La a 
EN u Hrn Ba Ka DE A DS SA DE 


Fo rs ToraTmafaorösraror 


3 Iteration 
44 SuM 
05 05 
rı SER 
i& A" 
42 STD 
ür 07 
43 REL 
oöi 01 
S5 + 

2 RÜL 
00 cu 
75 - 
3 RÜL 
[9 09 
55 + 
Ül 1 
ss = 
42 STO 
ü4 04 
37 1SZ 
03 09 
ss 18 
ie A" 
85 + 
43 RCL 
00 00 
15 - 
0i i 
935 = 
42 STD 
10 10 
71 SBR 
16 AR" 
a = 
01 1 
95 = 
42 STO 
07 07 
73 RC# 
7 07 
22 INV 
64 PI* 
10 10 
01 1 
22 IM 
44 SUM 
10 10 
43 RCL 
00 00 
85 + 
ci 1 
35 = 
E2 INV 
44 SUM 
7? 07 
43 RCL 
10 10 
53 - 
43 RÜL 
00 00 
95 = 
42 STO 


3.2 Inverse Iteration 69 


280 04 04 327 44 SUM 374 95 = 422 32 KT 
281 43 RCL 328 73 Rl# 375 42 STO 423 76 LBL 
282 10 10 323 07 07 376 09 09 424 53 <c 
283 85 + 330 22 INN 377 73 RC 425 53 « 
284 0 1 331 64 PD* se 10 10 426 73 Rl* 
285 95 = s3z 10 10 3793 50 IxI 427 10 10 
286 42 STO 333 43 RÜL se0 42 STO 428 50 IxI 
287 05 05 334 00 00 ss1 06 06 42939 75 - 
285 43 RCL 335 85 + 382 32 X1T 430 73 RC# 
289 00 00 336 01 i 383 76 LBL 431 07 07 
290 75 - 337 95 = 334 45 Yx 432 50 IxI 
291 01 1 338 22 INV 385 01 1 433 65 x 
292 95 = 339 44 5UM 386 44 SUM 434 43 RCL 
293 & STD 340 07 07 38?’ 10 10 435 06 06 
294 09 09 341 01 1 383 73 RC* 436 54 > 
295 7& LBEL 342 22 IN 339 10 10 437 50 IxI 
226 43 RÜL 343 44 SUM 390 50 Ix1 438 77 GE 
297 43 RCL 344 10 io 391 22 INV 439 54 9) 
238 00 00 345 43 RCL 392 77 GE 440 01 1 
239 75 - 346 10 10 393 52 EE 441 44 SUM 
300 43 RÜL 347 75 - 394 67 EU 442 10 10 
301 09 09 348 43 RCL 395 52 EE 443 44 SUM 
s02 95 = 349 00 00 396 42 STO 444 07 07 
303 42 STO 350 95 = 397’ 06 06 445 97 DS2Z 
304 08 08 351 42 STO 398 32 81T 446 09 09 
305 76 LBL 352 04 04 393 7& LBL 447 53 
306 44 SUM 353 71 SER 400 52 EE 448 61 GTO 
3507 73 RC& 354 16 A! 401 97 182 449 55 + 
308 04 04 355 85 + 402 09 09 450 76 LBL 
309 65 x 356 43 RÜL 403 45 YX 451 54 > 
310 73 Rl# 357 00 00 404 43 RÜL 452 97 DSZ 
311 05 0 358 75 - 405 00 00 453 02 02 
SIE 95 = 3593 01 1 406 42 STD 454 13 C 
313 22 IN 360 95 = 407 09 09 455 76 LBL 
si4 74 SM* 361 d2 STO 403 71 SBR 456 55 + 
315 10 10 3e2 05 05 409 16 A" 457° 25 CLR 
316 423 RCL 363 97 DS? 410 42 STO 458 91 R-S 
317 00 00 364 09 09 411 10 10 459 76 LBL 
sie 22 IN 365 43 RCL 4i2 85 + 460 16 A’! 
319 44 SUM 366 71 SBR 413 02 3 461 53 
s20 04 04 367 16 A! 4i4 65 x 462 43 RÜL 
3221 0 ı1 368 42 STO 415 43 RCL 463 00 00 
322 ZZ INV 369 10 10 416 00 00 464 33 x? 
323 44 SUM 370 43 RCL 41795 = 465 85 + 
224 05 05 371 00 00 415 42 STO 466 43 RCL 
325 97 DSZ 372 75 - 419 07 07 467 01 01 
32e 08 08 373 01 1 420 43 RCL 468 54 ) 
421 03 03 469 92 RTN 
Teil 3 
JO0 786 LBL nDü8 43 RLL 016 43 RÜL ü24 42 STO 
001 14 D 00% 01 Ui 017 00 00 025 05 05 
D02 43 REEL niü 95 = 018 95 = 026 76 LBL 
003 00 00 Oil 42 5T0 019 42 STD Ü27 65 
O04 zZ STD oiz2 10 10 020 r 7 0zE 73 Rlr 
305 03 093 113 85 + 021 00 0 029 10 10 
OlE 33 88 nid 02 2 022 42 STD 03033 88 
007 85 + D1S 65 8 023 0808 031 44 SUM 
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44 54 FI* (SEO ZE IMW 022 921 Re 
45 10 ı1n os? 42 FREI JE TE LEL 
code [1 1 os2 0505 rl 35 + 
O4 44 SM 053 43 FÜEL crl ’3 Rl# 
ads 10 An u un a nn orz jo 10 
242 44 SuM el 42 STD 072 1 RS 
so ur 209 08 74 ol 1 
cS1l 937 a 22 IMW O7S 44 5SLUM 
USE In b:| 2 44 51M TE 10 10 
053 65 rn U Eu üYT 27 DSE 
054423 + RCL DI = u 1 Do u 
055 08 SS 05 ao 0m + 
sn 21 R#S 


3.3 Der LR-Algorithmus 


Der Algorithmus von Rutishauser zur Bestimmung der Eigenwerte der regulären n,n-Matrix A 
beruht darauf, die Faktoren der LR-Zerlegung von A in umgekehrter Reihenfolge zu multipli- 
zieren und dieses Vorgehen zu wiederholen: 

A =: Ar = L; x R, 

R,'Lı=:Aa =L3'R 


RL; =:Ayyı = Liss  Risn 


Existiert die LR-Zerlegung jeder Matrix A; und sind alle Eigenwerte von A von verschiedenem 
Betrag, so konvergieren die L; gegen die Einheitsmatrix E und die R; gegen eine obere 
Dreiecksmatrix R, in deren Hauptdiagonalen die Eigenwerte der Matrix A stehen. 


Das Programm führt maximal eine vorzugebende Anzahl N von Iterationsschritten durch oder 
bricht vorher ab, falls 


IL; -EI, se lIAll, 


ist, wobei e eine vorzugebende Toleranzschranke ist. Vor Einlesen der Magnetkarten ist die 
Speicherbereichsverteilung mittels der Tastenfolge 3 2nd Op 17 auf 720 Programmspeicher- 


stellen zu ändern. Das Prograrnm bearbeitet Matrizen bis zur Ordnung n = 3. 


Programminstruktionen 


Toms Terme run] ame 


3 2nd 
Op 
17 719.29 
CLR 0 


Änderung der Speicherbereichsverteilung 


3.3 Der LR-Algorithmus 71 


Magnetkarten einlesen (Block 1, 2, 3) 
Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k = 15 


Eingabe der Matrix A spaltenweise 


6 | Ende der Koeffizienteneingabe 
Eingabe von n,k, N und e 


Anzeige der Eigenwerte A; 


Registerinhalte 


Roo: ---, Rıa: Programmzeiger 
Rıs, un. Rn2+1a: 411, -.-, Ann 
Rn2 +15, ---, Rntn-ı)+14: Zwischenergebnisse 


Bemerkung 


Ist eine LR-Zerlegung von A, nicht möglich, so hält das Programm und der Rechner zeigt 
dies durch eine blinkende Anzeige an. Nach Drücken der Tasten CLR und © erfolgt dann 
die Anzeige der Diagonalelemente von R;_;. 


Beispiel 


Gesucht sind in höchstens N = 5 Schritten mit der Toleranz e = 0.001 die Eigenwerte 
62 -1 

der Matrix A=|2 6 -3 
00 2 


72 3 Iteration 


Änderung der Speicherbereichsverteilung 


Magnetkarten einlesen (Block 1, 2, 3) 
Programmbeginn 
Eingabe von: k 


—_ 


441 
@21 
431 
412 
22 
432 
913 
923 
933 
Ende der Koeffizienteneingabe 


5 
6 
2 
0 
2 
6 
0 
1 
3 
2 


—_ 


2; 
OO OO NOT RWN- DW O 


Eingabe von: 
k 
N 
€ 


Anzeige von: X4 71.878787879 


4.121212121 
2 


-_ 
100 


Die exakten Eigenwerte sind A, =8, A2=4, A3=2 


Programm 3.3 Der LR-Algorithmus 
3 2nd Op 17 


oo0 76 LBL dizE 00 00 024 42 STD (32 91 RS 
oo 11 A 013 01 1 025 00 00 037 42 570 
002 91 R-S 014 44 SUM 026 42 STD 038 2 2 
005 42 STO 015 00 00 027 09 09 053 76 LEL 
004 00 00 üle 44 ZUM 025 91 R£S ü4lü ZZ IMY 

005 01 1 oir 01 01 029 42 ST0 041 43 RÜL 
006 42 STD ülS 61 GTO c3ca 0 01 odz2 00 00 
007 01 01 013 00 00 ü3l 42 ST0 042 42 5STD 
DS 43 RCL 020 03 05 ü32 02 03 ü44 05 05 
009 01 01 021 7& LBL ü22 91 R#S 045 Te LEL 
ll 91 RS ü22 12 68 034 42 5T0O ü4s 23 LH: 
011 72 ST# es 91 RS 1353 2 08 ü4?° 73 Rl# 


3.3 Der LR-Algorithmus 73 


O4 03 03 11lü Z2 IM ir2 73 Rl$ 
a493 50 I=I i1ll 64 FPI* ir2 04 04 
cSo 44 SUM liz 03 02 174 95 = 

Sl 11 11 112 01 1 175 22 INWY 
os2 01 1 114 44 SUM lfe& 74 SHE 
053 44 SUM 115 03 03 177 05 05 
0S4 02 03 ll& #7 ISE 175 43 RCL 
055 47 ISE ilF 04 04 179 00 00 
056 08 08 118 25 CLE 150 44 Sum 
(ST 23 LH# 119 43 RÜL 181 03 07 
055 43 RÜL 120 01 01 iS2E 01 1 

c59 11 11 121 85 + 153 44 Sum 
sc ” GE lzEE = RÜL 184 04 04 
lei 00 00 z2 00 00 135 97 152 
TS Se 124 95 = 15E 09 09 
cs: 00 0 125 42 STD 157 34 T= 

0e&4 42 5STD 126 04 04 1E8 43 RCL 
065 11 11 lET 85 + 182 10 10 
Üße 61 GTO 128 01 1 120 65 

ler 00 00 12399 = 191 43 RCL 
ü&S 7373 120 42 5T0 192 00 00 
064 32 81T 131 05 05 143 75 - 

070 00 0 122 42 STO 124 01 

Orl 42 STD 1322 0& 06 1945 95 

orz 11 11 1 R 136 82 

U7S3 97 1DSZ 13 147 44 

Dvd 092 09 1 19503 

075 ZZ IM 1 129 43 

Ür6 32 84T 1: z00 10 

077 42 STO 13 z0j 22 

ors 11 11 i zu: 44 

073 7& LBL 203 04 

080 24 CE 204 01 


jun) 
[ee] 
— 
Fu 
24 
a 
[er] 
Fr 


0853 95 = 2 

086 42 5TO 2 

087 04 04 149 32 84 273 

DEE 3 REL 150 3 RCL zr4 

os 01 0 151 00 00 z27E 

030 42 STO 15275 - = er6 

041 05 05 153 = RÜL 1 Eerr 

[32 89 + 154 07 07 = TE 

a3 01 15595 = STD 279 95 

094 99 = 156 42 STD lE ze0 22 

195 42 STO 157 08 08 Eu Zu] 44 SUM 
096 03 02 1Se 01 1 GE 282 03 03 
037 00 0 159 42 5T0 LEL } 1 
098 32 Kr lea 10 10 1.2 SLUM 
199 73 Rl# 1E1 76 LEL RÜL 

100 05 0 162 33 M8 OO 


102 87 EQ i64 10 10 


un ee HE SEEN um EEE Ener ES Au Ei Eu ne 
N CN n DE u ON u NE A) Kan ER ER Bin DE 5 
= 
ug] 
Fr 


[a Ed Du DE N 
2 Kur ur ı Ba Di 
NE MT (En mn 
- 

oO 
ea 
[77 
EN 
- 

x 


- 
c 
2; 

je Eu 5 Du ED EN LEN E 23 


107 ÖLR 169 73 RC# 
108 Ri Ir0 03 03 STD 
103 os 171 65 x Rc# 


ZH 
roao 


SZile+i 


oO 


rn 


rm“ 
EKD Qül4dle 


m DAY — 
a u WE Ws 
min 


TEL EN LEN 
en 
f 


a [rl 
m 
Dur Bu Bu 
17 


GHl+CcoMmo{oli- 
"MO 


Tl Fa 


Ar] 


OO 


233343 RCL 
334 42 RCL 
3535 08 08 
336 42 STD 
ss? 09 09 
s3e 76 LBL 
2339 44 SUM 
De} 


= e) 

+01° [a 

41l r3 Rls 
411 04 04 


[dr 
a u} 


m 
Dur Eu 2) 


D 
04 
L 
] 


d’co 00 00 
4ra =E IHM 


3 Iteration 


IE 


FEL 
[r] 


con 


Nr a 


3.4 Iteration in einer Variablen 


s44 43 STO 525  F2 Rl# | 95 
s45 04 04 SSE Üd e 95 = 42 
© 43 RÜL 587 42 STD 02 
cc sss y a2 SS 
- ssa m 1 F& LEL 43 
1 Sao 44 SUM ES + cc 
= Sal 02 02 43 RC 33 
STO S32 44 SuM 35 

s 

g. 


LEW} 
= 
] 
ım 


42 RFÜL rÜüe 
0 00 "Or 


= 
“maoqgi 


Co 
cı 
2 


42 RÜL s3 
cf ai Ir 
EI + 02 
od eo 54 
us = ed 
42 STD eo4 12 
03 03 Ben 42 
EEE + be) ESE Cd 
43 FÜL E ST 7 
co cn + BEE 04 
22 ME 4 eb u $ 
92 549 E20 01 
3 4? [1] 51 42 
\ E51 IE 053 
? E52 592 76 
F 553 E94 24 
r - Ed BIS 73 
r t de 5 855 E9E 03 
ar = 02 03 &SE Bar al 
Sr ° STO ES + Bar E95 01 
STE : 08 43 FÜL ESE E99 55 
ST’ TE LEL do 00 [te] on 3 
Sr 53 c 33 x8 ee ur 00 
Ss72 43 REL 95 = Bl "Oz 
seo 08 08 42 STD Ex "ü2 
sel 42 STO Ö4 04 4 
5 [13 3 42 
= & L 
5 4 


[Eu 


=] 
n 


=) 

— 

L 
"Leo Good Rica a 
— Peoalichgini 


3.4 Iteration in einer Variablen 


Gegeben sei eine kontrahierende Abbildung f, die der Lipschitzbedingung 
id -flyI<L-Ix<yl 


mit O<L<1 für x,y€[a, b] genügt. Das Programm liefert den Fixpunkt s = f (s) als 
Grenzwert der Folge 


xt); i=0,1,... 


12 


m+ 
vr 
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für jeden Startwert x, € [a, b]. Es bricht ab, falls die Höchstzahl N von Iterationen durch- 
geführt wurde oder falls für vorzugebendes e >O gilt: 
1-L 

>: 


Iui-xi-ılse' 


Die Abbildung f wird als Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift f(x) ist in Klammern einzuschließen. 

2. Für x ist RCL 00 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


1| Magnetkarte einlesen (Block 1) 


2| Eingabe der Funktionsvorschrift 


3| Programmbeginn 


4| Eingabe von xo,L, e, N 


5| Ergebnisanzeige 


Registerinhalte 


Roo, -.., Ros: Programmzeiger 


Beispiel 
Gesucht ist s= f(s) mit f(x) = = im Intervall [0,2], e= 10" und N = 10. 


Dabei sei L = 0.5, denn esgilt |f’(x)]| = | Zs | <05 


3.4 Iteration in einer Variablen 77 


Umschalten auf Bogenmaß 


Magnetkarte einlesen (Block 1) 


Eingabe von f(x) 


Programmbeginn 
Eingabe von: xo . 0.1 
0.5 
1 00 
1-10 
1-10 


.0000000001 
N 


Anzeige von: s .4501835576 


ao FE LEL 43 FÜL 6 LEL 42 43 REL 
ac dl A 2 12 5 43 04 04 
az 21 Rrs EI = r 44 42 STO 
al2 42 STD Ihe} D8 RR 245 dl an 
ao 00 a1 i 4 Idee 37 IS2 
als 31 R-S 75 - co4 04 247 02 05 
OO 42 STD +7 2 - 42 13 5 

Lu? 1 [2 43 RÜL 043 75 LEL 
Inn F- "5 4 oo an so 14 UL 

[alu b=) D aa + 33 = 21 42 FEL 
oll IE 43 FÜ so I=I 22 040 04 
cl 2 0: = 39 22 IHW 03221 Res 
c1E 33 = a4 To GE a4 FE LEL 


Kririe: 041 14 20 NS 22 BE 
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3.5 Steffensen-Iteration 


Gegeben sei eine in [a, b] kontrahierende Abbildung f. Dann konvergiert die Folge 
(f(x;) ei x;) 


EN re RR i=0,1,... 


KH KT 
gegen den Fixpunkt s=f(s) von f. Vorzugeben ist eine Toleranz n >O und die Höchst- 
zahl N von Iterationen. Das Programm bricht bereits nach weniger als N Schritten ab, wenn 


&Krı-l<n 
ist. Die Abbildung f wird als Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 
1. Die Funktionsvorschrift f(x) ist in Klammern einzuschließen. 


2. Für x ist RCL 00 zu setzen. 
3. Die Taste = darf nicht verwendet werden. 


4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift f(x} 


Programmbeginn 
Eingabe von x9, n, N 


Ergebnisanzeige 


Registerinhalte 


Roo, :.., Ros: Programmzeiger 


Beispiel 


Gesucht ist der Fixpunkt s=f(s) von f(x) = — im Intervall [0,3] mit n= 107°, 


N = 10 und dem Startwert xo = 0.1. 


3.5 Steffensen-Iteration 


Umschalten auf Bogenmaß 


Magnetkarte einlesen (Block 1) 


Programmbeginn 
Eingabe von: xXo 
n 


Anzeige von: s 


BIN IN! 
co 
DE 
DE 
aog 
JS 
IE 
007 
In] F=} 
cc 
10 
cl 
ulz 
013 
O1d4 
015 
OjE 
17 
[18 


Eingabe der Funktionsvorschrift f(x) 


4 
[: 


h | 30 

41 
04 
[NE 96% 
44 
45 
115 72 
IC 24 
IE: 32} 
243 
cS5ü 
021 
02E 
SE 
054 
0532 
056 


n 


RT 
20 


—Z 
Dr a ı BE N 


DU DE Dei Se 


n 
TI 


SF 
us8 
{137 
InT-A] 
Oel 
übe 


79 


1 


2 IHW 
ZUM 


jun] 


= REOL 


ci 


I=I 


IN" 
GE 
L 


" DSE 
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3.6 Das Newton-Verfahren 


Ist s einfache Nullstelle der differenzierbaren Funktion g, so existiert eine Umgebung von s, 
so daß die Folge 


si) 


ar) i=0,1,... 


Kr % 


für jeden Startwert xo aus dieser Umgebung gegen s konvergiert. Vorzugeben ist eine 
Toleranz n>O und die Höchstzahl N von Iterationsschritten. Das Programm bricht bereits 
nach weniger als N Schritten ab, wenn 

glx;) 
g’ (x) 


<n 


ist. Die Funktion g/g’ wird als Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift g(x) /g’(x) ist in Klammern einzuschließen. 
2. Für x ist RCL 00 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Ten Tee [Tone] an 


1 | Magnetkarte einlesen (Block 1) 

2 | Eingabe der Funktionsvorschrift g(x) /g’(x) 

3 |Programmbeginn 

4 | Eingabe von x9,n, N Xo 
n 
N 

5 | Ergebnisanzeige 


Registerinhalte 


Roo, ---, Ro2: Programmzeiger 


3.6 Das Newton-Verfahren 81 


Beispiel 
Gesucht ist eine Nullstelle von g(x) =x-cosx mit N=5, x, =0.5 und n= 10°'9, Es ist 


ax) _ X Zc0sX 
g’(x} 1+sinx 


Umschalten auf Bogenmaß 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift g(x) /g’(x) 


Programmbeginn 


Eingabe von: xo 
n 


N 
Anzeige von: s .7390851332 
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Programm 3.6 Das Newton-Verfahren 


aan Te ca cl 0 DT u In] l2 E 

ooi 11 ill Fa LEL 42 FEL TE LBL 
ans 51 ci 12 EB IE 02 13 5 

[u Tu RC Be: Be a1: 17 sa I=I 43 FÜL 
ang Od 012 Z2 InW Oo 00 
cos 21 o14 FT GE 1 Res 
Ode 32 als 02 02 12 5 "5 LEL 
Sur al Die 22 IMW ar ISEZ I ME 

[u Tu Po 9 217 44 ZUM ci al 


3.7 Regula falsi 


Ist s Nullstelle der Funktion g und liegen die Startwerte x, und x} genügend nah bei s, 
so konvergiert die Folge 


X %-1 


Kr Xn si; i=1,2,... 


glx;) -glx;_ı 
gegen s. Die Toleranz n >O und die Maximalzahl N von Iterationen sind vorzugeben. 
Das Programm bricht ab, wenn 
XKT%-1 
ee 
PER T EB 


ist oder N Iterationsschritte durchgeführt wurden. Die Funktion g wird als Unterprogramm 
eingegeben. Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift g(x) ist in Klammern einzuschließen. 

2. Für x ist RCL 00 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Te} 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift g(x) 


1 
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Term _ [ame Tr ae — 


Programmbeginn 


Eingabe von x%g, %1,7, N 


Ergebnisanzeige 


Registerinhalte 


Roo, --., Roa: Programmzeiger 


Beispiel 


Gesucht ist eine Nullstelle von g(x) =x -cosx mit xo = 0.4, x7=0.5, n = 10" 
und N = 10. 


Umschalten auf Bogenmaß 


Magnetkarte einlesen (Block 1) 


Eingabe der Funktionsvorschrift g(x) 


Programmbeginn 
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0.4 


0.4 
0.5 
1.00 
0 00 


Eingabe von: 


0 


Anzeige von: s .7390851332 


cool Fe LEL in 33 

ocoı 11 A U; +42 
als 21 Ra u cd 
Doz 42 STD ü Er 
IE C u 1] IE 
os 31 R-z iM 12 
Ale 42 5TO In) TE 
or 01 01 1 12 
oüE 21 Rs o 43 
003 32 MET 1) [1 
oil 91 5 [I £ EB 
il 42 0 ü 22 FE 

DU Be u Ü sc 77 GE 24 
012 FE LEL iu) 0S1 12.0 42 
14 12 8 [N GS2 43 RÜL co 
uüls 53 1 u IS 01 21 FE ! 
oile 53 X A] 54 55 + 32 BE 
ul? 42 FÜL ü 057 4: FEL 

als 01 o1 USE 0 03 [12 


3.8 Das Horner-Schema 


Zur Auswertung eines Polynoms 
pol (A) a0" ta," +... +an 
an einer Stelle X =, schreibt man zweckmäßigerweise 
p:= polo) = (...(a0Ao +24)  Antaa) An... +tamn_ı) Ant 


und arbeitet die Klammern von innen nach außen ab. Das Programm hat so nur Additionen 
und Multiplikationen mit dem festen Faktor A = Xu durchzuführen. Eingegeben werden nur 
die Koeffizienten a9, ..., an von pol {A). 


Beachte: ao ist der Koeffizient von A"! 
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Magnetkarte einiesen (Block 1) 


f) 
Programmbeginn DER ? 


Eingabe von n; a9, 41, ..-, An 


Ende der Koeffizienteneingabe Iren { 
Eingabe von Ay 
Ergebnisanzeige 


eine Jah 


Registerinhalte 

Roo, :--, Ros: Programmzeiger 
Rose: LEST Rn+6: 90, 41, .--, An 
Beispiel 


Zu berechnen ist der Funktionswert des Polynoms 
pol(A)=-M+3MR +41-2 
an den Stellen X, =0.5 und X, = 1.5 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von: n 
a0 
a1 
a2 
a3 
Ende der Koeffizienteneingabe 


BR2AOD-O -— 


Eingabe von: %o 

Anzeige von: p, = pol (A,) 
anderes Argument 

Eingabe von: A, 


Anzeige von: p, = pol (X,) 
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Programm 3.8 Das Horner-Schema 


InInIE) ve LEL 
oil Al A 
PIE a1 RES Von? 

aos 42 5STD or 
Du Fe a un ZUM 27 
ads ad u 05709 
a0 42 5TO Zum 
U u Pe u 2) 
[ETF ur E 

20a 42 5T0O 

iu 05.05 

ll TE LEL 


m 


3.9 Das erweiterte Horner-Schema 


Pau sale arreHll 


I 

42 sro 
24 ar! 

43 FÜL 

oc 


[45 
Ode 
Nor 


N 
een 


[N 
% 


Das erweiterte Horner-Schema liefert außer dem Funktionswert p des Polynoms 


po) = ta +... +2, 


01° 1 
44 SM 


(siehe 3.8 „Das Horner-Schema’’) auch den Wert der Ableitung q = pol’ (X) an der Stelle 


A = ‘og. Eingegeben werden nur die Koeffizienten a9, ..., an von pol (A). 


Beachte: ao ist der Koeffizient von X"! 


Programminstruktionen 


7 Tran Tone [rem nie | 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von n; ag, &4, ..-, An 


Ende der Koeffizienteneingabe 
Eingabe von Xo 
Ergebnisanzeige 


Registerinhalte 


Roo: -.., Ros: Programmzeiger 


Roz, ... 


P Ru+7: a0, 91, --- 


ı An 
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Beispiel 

Gesucht sind Funktionswerte und 1. Ableitungen des Polynoms 
pol (A) =-M+3R +4\-2 

an den Stellen Ag =-1.5 und A,=2.5 


Anmerkungen Eingabe 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von: n 


a0 
a] 
a 
a3 
Ende der Koeffizienteneingabe 


PB RROND-O — 


Eingabe von: Ay 


2.125 
R/S -11.75 


Anzeige von: p,= pol (Ay) 
Go = pol’ (Au) 
anderes Argument 


Eingabe von: A; 


Anzeige von: p,=pol(A}) 
9,” pol’ (A) 


cool Te LEL Di 1 03 02 D55 
Sol dl A dd Sum In USE 
Doz 21 R#s DE 06 42 STD 57 
Goa 42 STO 44 = > 15 LE Jg 
oc ac 02 08 os 43 RÜL 57 
OS u E1 5TO 41 ? DEN 
LNCHE, S Sa NE üg4z Del 
Sc } 4 IE 
005 12 \ 44 (IE: 
Dan be) 3 0145 De 
O1u 4: ode DES 
coli IR IE 47 Ei 
OLE 42 RÜL 4 ? 
u 9 Br Tu WIN! Ic 55) 1 
o14 42 STD JSc a 
CS oil 01 ul ü 
IB E= un TE 1 
fr 42 ST0 053 2 

4 FR 
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3.10 Einfache Nullstellen von Polynomen 


Einfache Nullstellen s eines Polynoms 
pol (A) =" ta" "+... +2, 


bestimmt man durch das Iterationsverfahren von Newton (siehe 3.6 „Das Newton-Verfahren”): 


pol (X,) . 
Auı= dr Bor)” i=0,1,... 
Dabei ist Au ein geeigneter Startwert. Der Quotient 
IR pol (A,} 
IT Bor (A) 


wird im erweiterten Horner-Schema bestimmt (siehe 3.9 „Das erweiterte Horner-Schema”). 
Eingegeben werden neben den Koeffizienten a9, ..., an und dem Startwert X, eine 
Toleranz e >O und die Höchstzahl N von Iterationen. Das Programm bricht nach weniger 
als N Schritten ab, wenn 


lal<e is. 


Beachte: a, ist der Koeffizient von X"! 


Programminstruktionen 


Teen Ferse Tree] ame 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von n, €, a0, :.-, An 


Ende der Koeffizienteneingabe 
Eingabe von X, und N 


Ergebnisanzeige 


Registerinhalte 


Roo. ..., Roz: Programmzeiger 
Ros. Per Rn+g: ag, 41, -.+, An 
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Beispiel 
Gesucht sind die Nullstellen s; , s2, s3 von 
pol(A)=-N+3MR +41 -2 
und den Startwerten Ao1=- 1, Aoa = 1, Aoa=5 und e= 10-10 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von: n 
€ 


a0 
a7 
92 
a3 

Ende der Koeffizienteneingabe 


I 
N 


Eingabe von: 


-1.292401585 
-1.292401585 
1 


Anzeige von: sı 


anderer Startwert 


Ao2 
N 


Eingabe von: 


‚3972950693 
‚3972950693 
5 


Anzeige von: 5 


anderer Startwert 
Eingabe von: Ag a 
N 


3.895106516 


Anzeige von: s3 


ct 


Beben m 


Far BT fa 
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“1 2 zz IMY 
+ TG " 

DE ds 

TE LEL 

13: IE 

42 RÜL 2 r 
TUE EIN; - r 
42 STD DEE 
ci 01 a7 DIESE 
un) lb DE 
42 STO 13 0 

202 TE LEL 
39 14 2 

42 =TD 42 Fü 

os 05 N 

42 FÜL 1 RS 


3.11 Das Verfahren von Bairstow 


Hat das reelle Polynom 
pol) = ta +... ran 


die komplexe Nullstelle u, so ist auch die zu u konjugiert komplexe Zahl & Nullstelle von 
pol (A) und (A-u) (A-Z) ist ein reelles quadratisches Polynom A -uA -v, das sich nach 
dem Euklidischen Algorithmus von pol (A) abspalten läßt. Sind uo und vg Näherungen 

für u und v in ® -uX-v, so liefert das Programm verbesserte Näherungen u; = Uo + Au 
und v7= vo + Av. 


Beachte: a9 ist der Koeffizient von A"! 


Programminstruktionen 


Verfahren Anzeige 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von n, ug, Vo, &0, 1, --- 


Ende der Koeffizienteneingabe 


Anzeige von u} und v; 
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Registerinhalte 
Roo: ---, Rıo: Programmzeiger 
R;1. ee Ra+n:: 80, --+, An 


Rn+12, -.-, Ran+ 1a: Zwischenergebnisse 


Bemerkung 


Das Verfahren ist nicht auf Paare konjugiert komplexer Nullstellen beschränkt. 


Beispiel 
Die Näherung des quadratischen Faktors 


KR -WA-vo=R-1.8X-(-2.3) 
von 
pol (A) = M-R-6%+14X -12 


soll verbessert werden. 


Anmerkungen Eingabe Anzeige 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von: n 
Uo 
Yo 
a0 
a; 
a, 
a3 
a4 
Ende der Koeffizienteneingabe 
Anzeige von: u; 1.941934625 
v4 -1.983128031 


GPWOND-O0D0Rn- - 


Esist u=2 und v=-2. 
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Programm 3.11 Das Verfahren von Bairstow 


Tr 
A) 


Eu 


m 
Da RL 


) 


ıTı 


aac Te LEL 108 Rü# E33 TE Als 
oo I A 103 nk} DE u re u Ir 
Ü 1 Res 110 = BES TS - 
co3 42 5TO ST EEO TS Rü 
oo 10 1u - = 

4 £ Pe 

4 e 


n 
4 
3 
2 
1 


111 
[ EZ 
5 113 3 
u [a 114 | E S 
[N In] co 3 119 : Be 
[N = es 1 11 ii 05 
Ixnke] 2 117 Dr Kr 
ac E IE =TDO 118 73 > Fl# 
c1l : ute] 11? 4 06 06 
I1E N Re 120 De =: 
[1: ES = E YET 
214 43 FL “ro 07 
DT un) 05 re 7 
85 + Do 
cl F Rüs 20 0 
O4 IR] 24 be & 
35 = = F 
42 - :T# £ 
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(u He) 
a a u On ı BR Cu 
Dr 


an 
{Ho l4oal po + 
oO 


EN EN EL lan 


[7 I Ü a A 

Iuln] 1 5 r 
TE STE } 42 FELL 
’ dr 1) ne! 

: REL a = 
ur zE InW 

+ 20 44 Sup 
1 ai 01 ül 

b = 32 HE Dr 
= 5STD 33 T2 RC# 
5 05 a4 030693 

+ 
1 1 3 72 Re 
SS = sr 04 04 

= 5STD \ A 
hr Er] S v2 Rl# 

Bin 0 u Ede 

Ss + = 5T# EM 
1 1 05 ra Rl#$ 
odlo 023 EL E08 
ü4dl 35 = 5 54 
de 42 STO be 

uf: a 03 43 


244 
145 
Ud46 


Din} 
- 
oO 


101 
108 
103 
104 
105 
Lüe 


Hl + 
3 


= 


un 
nn 


E 

n 
S 
= 
1 
ae? 
3, 
re} 
[4 
S 


+ To 


107 85 + 2 
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3.12 Das Bernoulli-Verfahren 


Ist pol (A) =A”+a,X” "+... +2, ein normiertes Polynom (a, = 1) vom Grad n, so ist 
pol (A) das charakteristische Polynom der Matrix 


0 1 
0 
und es gilt 
0 
0 1 
an a2 a 
0 1 Yk Ye+1 
0 
_ mit 
0 
0 1 Yk+n-2 Yk+n-1 
an » + 092 a Yk+n-1 Yk+n 


n 
Ya Da a k=0,1,... 


i=1 


Wählt man einen geeigneten Startvektor [yJ, ---. Yn_ı F; etwa [O,....,O, 1j"; und besitzt 
pol (A) eine betragsgrößte Nullstelle A}, so konvergiert die Folge der Quotienten Yerı/Ya 
gegen A}. Das Programm benutzt den oben angegebenen Startvektor und bricht nach der 
vorzugebenden Anzahl N von Iterationsschritten ab. Es bearbeitet auch Polynome mit au # 1 
und eignet sich zur Bestimmung eines Startwertes für das Programm 3.11 „Einfache Null- 
stellen von Polynomen’. 


Beachte: a, ist der Koeffizient von A”! 


Programminstruktionen 


Teanen  Tenese [re] ame 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von N, a0, 4, ..., An 


Ende der Koeffizienteneingabe 


Anzeige von A, 
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Registerinhalte 


Roo:; ---, Rog: Programmzeiger 
Ro7: ---, Rn+7: 80, 81, ---, dn 
Ra+ßı «+, Ran+7: Ykr 3 Yk+n—1 


Beispiel 
Gesucht ist in N =5 Schritten die betragsgrößte Wurzel A, von 
pol (A) =-N+3R+41-2 


Magnetkarte einlesen (Block 1) 


Programmbeginn 


Eingabe von: N 
a0 
a7 
a 
a3 
Ende der Koeffizienteneingabe 
Anzeige von: X; 3.880829016 


Programm 3.12 Das Bernoulli-Verfahren 


Scan FE LEL In zuge) 
iA ada 0 
1 FR#3 sol 7 ; 
42 STD sl 0 
un ET) USE TE LEL 
oo u 0532 Br 
42 STO S4 0 h 
O1 01 055 4 ; 
“ 5 USE 0 DE [ned = 
5 - IS? 43 RÜL SLIM 
210 1 osz 01 0 [6 
all = 153 859 8 1 
LE STD uE ae SLIM 
01: 0: E + DE 
id LEL IE " 
715 Ts (IE = 
Üle LE =TD 
Ulr de: [3 
le a7: u 
J1% LE STD 
zn IT: DE 
ü REL 
AT 1 
[M] S5TO 


(de) 
1 


3.13 Das inverse Bernoulli-Verfahren 


43 FRI ic I: a2 75 - 
DE 06 1il 1: ; uB 1 
42 RÜL 112 1: b ] = 
ci 01 112 1 di 42 ST 
3 - 114 1 id2 u [nel 
ol 1 115 12 i42 L 
SE; „= te 13 

Bun 42 370 11? 1 

cd ds nos 115 1 

05 585 + il? 1 

ice 095 120 1 

dr 99 = 11 1 

10a de STO l 

FON ee 122 13 


3.13 Das inverse Bernoulli-Verfahren 


Ist An #0 betragsmäßig kleinste Nullstelle von 
poA)=oN ta "+... +an, 

so ist Un = 1/A,, die betragsmäßig größte Nullstelle von 
rez (x) = u” pol 4) Sotraut..t+anM”. 


Wendet man die in 3.12 „Das Bernoulli-Verfahren’’ geschilderte Methode auf rez (x) an, 
so konvergiert die Folge der Quotienten y,/y«,. gegen An. Die Anzahl N der vom Pro- 
gramm durchzuführenden Iterationen ist vorzugeben. 


Beachte: ao ist der Koeffizient von A"! 


Programminstruktionen 


Magnetkarte einlesen (Block 1) 
Prograrnmbeginn 
Eingabe von N, ao, 44, ..., an 


Ende der Koeffizienteneingabe 


Ergebnisanzeige 
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Registerinhalte 


Roo, :--, Roz: Programmzeiger 
Ros. eg Rn+7: 80, -+-, n 
Rn+s, ... Ran+7: Yeı ... Yk+n-1 


Beispiel 


Gesucht ist in N = 5 Schritten die betragskleinste Wurzel von 
pol (A =-N+3% +41 -2 


Magnetkarte einlesen (Block 1) 


Programmbeginn 
Eingabe von: N 
a0 
a1 
a2 
a3 
Ende der Koeffizienteneingabe 
Anzeige von: X, ‚398255814 


au Fe LEL o4E 2 08 
ccoı 11 A Ü: 2143 £ 

o02 31 Rs Ü; oso r3 

ca 42 STO nF 0Si 03 

NINE um 00 [; [Uwe 35 = 
cas 000 Ü; 053 44 ZLM 
ÜDE 4 5TO Ü% a4 O6 Us 
oor 01 Ü: 055 [I 

NR To u Io Ü% ISE 

02 42 5TD [N UST 
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ulE IuT-I8] 


Tı 


A 
42 5STD 
3 


1 
IE 
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ÜBER U | 
er 42 5T0O 
oe DE DE 
E23 TE LEL 
rc 42 5STD 
rl 72 Rl# 
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a4 +5- 11ld 42 =T0O 124 44 SuM 
22 IHMW ill 02 08 125 02 02 
+42 FEI ilE 535 + 12E& 44 SUM 
uT23 DE 113 1 1F7 uReH 02 
+2 FÜL 114 = 122 Ar DSE 
O1 115 ST 123 05 05 
vs - 118 13 120 42 FRÜL 
I Es 11? ! 121 42 FÜL 
5 = 113 1 u u 1 
42 510 117 1 5 
0505 120 1: 1 
SS + 121 1% 

a5 1: 1: 

5 = 123 1: 


3.14 De OD-Algorithmus für tridiagonale Matrizen 


Ist A eine tridiagonale Matrix, deren obere Nebendiagonale aus Einsen besteht 
q, 1 


e,9, +, 1 


0 


en-19n-1 en-ıt An 
und existiert die LR-Zerlegung von A, so ist 


1 g, 1 


en-ı 1 "m 


Berechnet man nach den Regeln 


eli+m + ee e el) + au 


k-1 
und 

en = ei) a 
die Folgen a „., und (el?) - ‚ so konvergieren die e, gegen Null und die q, gegen 


die Eigenwerte von A, falls die anwerte von A paarweise von verschiedenem Betrag sind. 
Das Programm führt eine vorzugebende Anzahl N von OD-Schritten durch; dann erfolgt die 


Ausgabe von Een BT So ; ne u a: Wird im Verlauf der Rechnung eines der q, zu 


Null, so hält das Pranam und der Rechner zeigt dies durch eine blinkende Anzeige an. 
Die Ausgabe der momentanen e, und q, erfolgt dann nach Drücken der Taste C. 
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Programminstruktionen 


Programmbeginn 
Eingabe von e;, ..., en-1; 4, --- 


4 | Ende der Koeffizienteneingabe 
Eingabe von n und N 


Ergebnisanzeige 


Registerinhalte 


Roo, ---, Ros: Programmzeiger 


Rose, ne Rn+7: @1,..., &n—1 
Ru+8, or Ran+7: Q, +++, An 
Beispiel 


Gesucht sind die Eigenwerte der Matrix 


110 ı100|/|1 10 
A=|1 2 1|=|1 1 0|/|0 1 1|=LR. 
023 02 1)[0071 


Esist e,=1, e=2, q,=0=q,=1. Es sollen 20 OD-Schritte durchgeführt werden. 
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Anmerkungen 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von: e; 
e 
4, 
% 
93 
Ende der Koeffizienteneingabe 


oO oOQRWON  - — 


Eingabe von: 


Anzeige von: ‚0000005666 
5.2804618 -22 

4.114906557 

1.745898729 

.1391941469 
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5 
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ci 01 44 07 Fr 4 15 
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105 TE LEL 12 i4i 
1C3 IA je i 142 
11u 25 CLR 1 i4: 
111 02 8 1 144 
11z BEIM i 145 
115 43 FÜL 15: 146 
114 DO u E UT 1 147 
115 rar 1 1423 
116 ol 1 1 i43 
117 = 1 

115 42 =TD 1: 


3.15 Der OD-Algorithmus für Polynome 


Beginnt man zu einem gegebenen Polynom 
pol (A) = a0" ta, +... +5; 3%0 


das QD-Schema mit der horizontalen Doppelzeile 


N 


und bestimmt beginnend in der rechten oberen Ecke die Schrägzeilen nach den Rhomben- 
regeln (siehe 3.14 ‚Der QD-Algorithmus für tridiagonale Matrizen’’), so hat die zu einer voll- 
ständigen Schrägzeile gehörige Tridiagonalmatrix das gegebene Polynom zum charakteristi- 
schen Polynom. Das Programm führt nach Bestimmung der ersten vollständigen Schrägzeile 
eine vorzugebende Anzahl N von OD-Schritten durch; dann erfolgt die Ausgabe von 

N, eg. Sind die Nullstellen von pol (A) paarweise von verschiedenem 
Betrag, so konvergieren die e, gegen Null und die q, gegen die Nullstellen von pol (A). 

Wird im Verlauf der Rechnung eines der q, zu Null, so hält das Programm und der Rechner 
zeigt dies durch eine blinkende Anzeige an. Die Ausgabe der momentanen e( ei) 


! | 1 7 En, 
a", es gq(’ erfolgt dann nach Drücken der Taste C. 


Beachte: a9 ist der Koeffizient von A"! 


Programminstruktionen 


Magnetkarte einlesen (Block 1,2) 
Programmbeginn 


Eingabe von a9, a4, ..., An 
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Team —_ _____ Tenee [rau] am 


Ende der Koeffizienteneingabe 


Eingabe von n und N 


Ergebnisanzeige 


Registerinhalte 

Roo; ---, Rıo: Programmzeiger 
Rı1, ---, Rn+1:: a9, 91, ---, @n 
Ra+12, +--, Ran+1o: E1,.--, En—-ı 
Ran+11r +, R3an+10: 1, ++, 


Beispiel 
Es sollen mit N=5 OD-Schritten Näherungen für die Nullstellen von 
pol (A)=-N+3M+4\-2 


bestimmt werden. 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 
Eingabe von: ao 

a1 

a2 

a3 


Ende der Koeffizienteneingabe 


OP POND- OD 


Eingabe von: n 
N 
Anzeige von: e, -.0063368755 
e) -.0003864977 
q, 3.899866489 
% -1.290529036 
Ga ‚3973859207 
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Programm 3.15 Der OD-Algorithmus für Polynome 
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3.15 Der OD-Algorithmus für Polynome 
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4 Interpolation und diskrete Approximation 


4.1 Lagrange-Interpolation 


Zu gegebenen Stützstellen xo, ..., Xn bilden die Lagrange-Polynome 


ge a en 2 mL 
ER: X)... Xi) Kirn) ee (RT nn) 


eine Basis im Vektorraum der Polynome bis zum Grad n. Es ist 
Lu)=8 | 1 m i=k 
. x = . = ur 
BUNS ae i#k 


Sind zu den Stützstellen xg, ..., Xn Stützwerte fo, ..., f„ gegeben, dann hat das Inter- 
polationspolynom durch die Knoten (x,, f;) die Form 


pol =fo lol t... ta intX), denn es ist 
pol (x.) = 2 to) = B fh ö-fe: 
i=0 i=0 


Das Programm liefert den Wert des Interpolationspolynoms an der Stelle x durch Auswerten 
von (*). 


Programminstruktionen 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von n, Xo, --», Xns for - +, fn 


n+2 


2nt2 
Ende der Koeffizienteneingabe 2n+2 
Eingabe von x 
Ergebnisanzeige pol (x) 
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Registerinhalte 

Roo. --., Ros: Programmzeiger 
Ro7, -.-, Rn+7: Xo4 +:-+ Xn 
Ru+8: es Ran+e: fo, ln 


Ran+9ı ---, Ran+g: Zwischenergebnisse 


Bemerkungen 


1. Die Schritte 4 bis 6 der Programminstruktionen lassen sich mit beliebig vielen 
Argumenten x wiederholen; 

2. Wird x =x; (also eine der Stützstellen) als Argument eingegeben, so hält das Programm 
und der Rechner zeigt dies durch eine blinkende Anzeige an. 


Beispiel 
E x 0 1 2 } 
Das durch die Tabelle ea. gegebene Interpolationspolynom soll an den 
i 
Stellen x =3 und x’ = - 1 ausgewertet werden. 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von: n 


1 
1 
0 
1 
2 
3 
4 
5 
6 
6 


. Ende der Koeffizienteneingabe 
Eingabe von: x 
Anzeige von: pol (x) 
anderes Argument 


Eingabe von: x 
Anzeige von: pol (x’) 
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Programm 4.1 Lagrange-Interpolation 


LEL ze + ee] 
A [u [u ei 
RS 25 = Sa 
STO 42 ST0 [13 
ac 02 113 39 
5 ur 7 +42 
=T0O 42 STD S 
1 O4 04 
TE LEL 
a4 TH 
42 FRÜL 
cd InIR] 
+ 
1 
STD 
ne 
1 
=T# 
0: 
STO 
.S 
LEL 
BE 1-8 
VE j n 
12 ! 5 
4 1 Bo E “Ti 
42 57T FE 110 am cm 
IE u 73 Rl# ill 85 + 
42 Fü 05 05 112 oı 4 
ao a S4 1l2 35 = 
25 + Er EU 114 42 570 5 
ci 41 10 E' 115 01 01 ZE 
35 = E4 FD* lie 35 + 197 
42 STO 02 02 llr 077 132 
Bi TE LEL lin 90 = 153 
: 10 E* 11# 42 57T 1E0 42 FÜL 
ci 1 lz0 02 08 1E1 04 0üg4 
44 SUM 43 R lEE 91 RrS 
Jg 05 05 DC 


RP, RX) Pisa, 


u, KIT 0 TFT, u 


X = XXX 


Dabei ist p, = f;. Das Programm bestimmt den Wert des Interpolationspolynoms an der 
Stelle x nach dem Schema von Neville. 
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Programminstruktionen 


1 | Magnetkarte einlesen (Block 1) 


2 | Programmbeginn 
3 | Eingabe von n, Xg, :.-, Xnı fo, ---, fn 


4 | Ende der Koeffizienteneingabe 
5 | Eingabe von x 


Registerinhalte 


Roo, ---, Roz: Programmzeiger 
Ros: ER Rn+g: X0, ++, Xn 
Rn+9, eg Ran+9: Po: Kine Pı 


Bemerkung 


Da die Konstanten fo, ..., fn „überschrieben’’ werden, läßt sich das Programm nicht zur 
Auswertung des Interpolationspolynoms an mehreren Stellen verwenden. 


Beispiel 
. xx|0o 1 2 ö 
Das durch die Tabelle : ae gegebene Interpolationspolynom soll an der 
i 
Stelle x = 3 ausgewertet werden. 
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Magnetkarte einlesen (Block 1) 


Anzeige 


Programmbeginn 


Eingabe von: n 


Ende der Koeffizienteneingabe 
Eingabe von: x 


a u 


Anzeige von: Po 1.2 


Programm 4.2 Das Schema von Neville 
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4.3 Entwickeln nach Tschebyscheff-Polynomen 


Die Tschebyscheff-Polynome T,„ (x) bestimmen sich rekursiv aus den Formeln 


Tolx)=1; Tl) =x 
Tn+1X) = 2x Talx) - Tn-ılX) . 


Da von allen normierten Polynomen vom Grad n > 1 das Polynom 2'""T, (x) in [-1, 1] 
die kleinste Tschebyscheff-Norm besitzt, ist es nützlich, ein gegebenes Polynom 


pol(x) = co +cıX +... +c„x" 
in Tschebyscheff-Polynomen zu entwickeln, also zu schreiben 
pol (x) =aoTolx) +a1 Tılx) +... tan Tn (x). 


Will man pol (x) nämlich durch ein Polynom vom Grad n-1 annähern, so läßt man das 
letzte Glied der Tschebyscheff-Entwicklung fort; der maximale Fehler ist dann wegen 
IT, = 1 auf [-1, 1] gerade |a„|. Das ist die nach der Tschebyscheff-Norm bestmögliche 
Approximation eines Polynoms vom Grad n durch ein Polynom vom Grad n -1. 

Das Programm entwickelt Polynome bis zum Grad n = 15 in Tschebyscheff-Polynomen 
und gibt dann die Koeffizienten a; aus. 

Um dieses Programm auf Magnetkarten zu speichern, müssen außer den Programmschritten 
auch Daten, nämlich die Koeffizienten der Tschebyscheff-Polynome, aufgezeichnet werden. 
Dazu geht man wie folgt vor: 


1. Änderung der Speicherbereichsverteilung auf 100 Datenspeicher mittels der 
Tastenfolge 10 2nd Op 17. 

2. Eingabe der Programmschritte (siehe Programmausdruck). 

3. Abspeichern der Konstanten in den jeweiligen Datenspeichern (siehe Ausdruck 
der Registerinhalte). 


4. Beschreiben der Magnetkarten (Block 1, 2, 3, 4). 


Programminstruktionen 


Änderung der Speicherbereichsverteilung 


Magnetkarten einlesen (Block 1, 2,3, 4) 
Programmbeginn 


Eingabe von Co, Cı, .--, Cn 
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Ende der Koeffizienteneingabe B 
Ergebnisanzeige 
R/S 
R/S 


Registerinhalte 
Roo; -.., Rog: Programmzeiger 
Rio, a Rn+10: Co, »..,En 


Rag, ..., Rgg: Koeffizienten der Tschebyscheff-Polynome 
(siehe Ausdruck der Registerinhalte) 


Beispiel 


Das Polynom pol (x) = 2 — 9x + 4x? + 16x? + 8x* soll nach Tschebyscheff-Polynomen 
entwickelt werden. 


Änderung der Speicherbereichsverteilung 


Magnetkarten einlesen (Block 1, 2,3, 4) 
Prograrnmbeginn 
Eingabe von: Co 

cı 

c2 

c3 


run. 0 oo %»20 


Ca 
Ende der Koeffizienteneingabe 


Anzeige von: &0 
a1 
a2 


a3 
94 


Die gesuchte Tschebyscheff-Entwicklung lautet also 
pol (x) =7'Tolx) +3-T,(x) +6 Ta(x) +4 Tal) + Tale). 
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Entwickeln nach Tschebyscheff-Polynomen 
10 2nd Op 17 
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SE Er 23 
Sr BE - Ein 
SE 53 Ep 
3 rÜü 32 
2) A 33 
#1 TE a4 
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55 FE -Eidal, 98 
203 Fr 1E204 ER 


4.4 Ökonomisieren eines Polynoms 


Die Tschebyscheff-Entwicklung eines Polynoms ist besonders nützlich, um seinen Grad zu 
ökonomisieren. Wegen ||T,|I, = ? in [-1, 1] ist der maximale Fehler, der durch Fortlassen 
des letzten Gliedes entsteht, höchstens |a„,|. Will man das Polynom 


pol(x)=cotcıxX +... tc„xX" 


durch ein Polynom möglichst niedrigen Grades approximieren und dabei höchstens den 
Fehler € >O begehen, so streicht man in der Tschebyscheff-Entwicklung 


pol (x) =a0To(x) +... +a,T, (x) 


solange das jeweils letzte Glied, bis 


n 
D, lal>e 
i=k 


ist. 
Das Approximationspolynom ist dann 


n k 
app (x) = pol (x) - > a,T;(x) = 2 b;x! 
i=k+1 i=1 
mit dem Grad k. Das Programm liefert die Koeffizienten b; von app (x) in der gewöhn- 
lichen Basis 1,x, ..., x*. Es ökonomisiert Polynome bis zum Grad n = 15. 
Um dieses Programm auf Magnetkarten zu speichern, müssen außer den Programmschritten 
auch Daten, nämlich die Koeffizienten der Tschebyscheff-Polynome, aufgezeichnet werden. 
Dazu geht man wie folgt vor: 


1. Änderung der Speicherbereichsverteilung auf 100 Datenspeicher mittels der 
Tastenfolge 10 2nd Op 17. 
. Eingabe der Programmschritte (siehe Programmausdruck). 
3. Abspeichern der Konstanten in den jeweiligen Datenspeichern (siehe Ausdruck 
der Registerinhalte). 
4. Beschreiben der Magnetkarten (Block 1, 2, 3, 4). 


[8 
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Programminstruktionen 


[Tom _______ [ense[ re] Aue 


1 | Änderung der Speicherbereichsverteilung 


Magnetkarten einlesen (Block 1, 2, 3, 4) 
Programmbeginn 


Eingabe von e&, Co, ..-, Cn 


Ende der Koeffizienteneingabe 
Anzeige von k = Grad app (x) 
Anzeige von bo, ..., bk 


Registerinhalte 


Roo, -.., Rog: Programmzeiger 

Rıo, ee Rn+10: Co: +++, En 

Ra7: &; 

Ras, ..., Rgg: Koeffizienten der Tschebyscheff-Polynome 
(siehe Ausdruck der Registerinhalte) 


Beispiel 


Das Polynom pol (x) = 1.571 x - 0.646 x?+ 0.08 x® soll durch ein Polynom app (x) 
niedrigeren Grades approximiert werden, so daß 


IIpol (x) = app (x)Il, < 0.01 


ist. 
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Änderung der Speicherbereichsverteilung 


Magnetkarten einlesen (Block 1, 2, 3, 4) 
Programmbeginn 
Eingabe von: € 

Co 

cı 

a 

c3 

Ca 

Cs 


0 
4 
0 
0 
1 
2 
3 
4 
5 
6 


Ende der Koeffizienteneingabe 
Anzeige von: k 


b, 


Also ist app(x) = 1.546 x - 0.546 x?. 


Ökonomisieren eines Polynoms 
10 2nd Op 17 
Programmteil 


i DT ut) 
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4.5 Methode der kleinsten Quadrate 


Bei der diskreten Approximation der Funktion f durch ein Polynom vom Grad n im 
Intervall [-1, 1] wird die Genauigkeit möglicherweise erhöht, wenn man die Anzahl der 
stützenden Punkte (x;, f}) auf m+ 1 mit m>n erhöht. Wählt man die Tschebyscheff- 
Polynome To (x), ...., Tn (x) als Approximationsfunktionen und sind die Stützstellen 
XO4 +++, Xm die Nullstellen des Tschebyscheff-Polynoms Tm+ (X), so ergibt sich wegen 
der Orthogonalität der Tschebyscheff-Polynome das folgende besonders einfache 
Gleichungssystem für die Koeffizienten a; der Approximation 

pol(x)=agTolx) +... tan Tn (x): 


f (xo) 
m+t1 a9 Tolxo). 5 . To (Xm) 


0 m+1||a, Ta&xo): . . TnXm} 
2 f(x) 


in der Tabelle unten sind diejenigen Paare (n, m) mit + gekennzeichnet, für die der Speicher- 
platz ausreicht. Wird der Rechner mittels der Tastenfolge 7 2nd Op 17 auf 70 Daten- 


speicher umgeschaltet, wird die Approximation auch für die mit o gekennzeichneten (n, m) 
berechnet. 


Magnetkarte einlesen ivlock 1, 2) 
Programmbeginn 


Eingabe von n und m 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k > 12 
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5 | Eingabe von X%og, .--, Xmı for ---, fm 


6 | Ende der Koeffizienteneingabe 
7 | Ergebnisanzeige 


Registerinhalte 


Roo, -.., Rıı: Programmzeiger 

Rx, .--, Re+n: &0, ---, an 

Re+n+1, :--, Rktnt+m+1: X0s +, Xm 

Resntm+2, +, Re+n+2m+2: fo, ---, fr 

Re+n+2m+3: "u Re+2n+3m+nm+3° Tol&o), ..., Tn(Xm} 


Bemerkungen 


1. Ist f(x) nicht auf [-1, 1], sondern auf [a, b] definiert, so ist f(x) nach f(t) mit 
ge 2x-a-b 
b-a 
zu transformieren. 
2. Die Nullstellen von Tn+1(x) sind 
2itim _ 
m+1 2: J 


x;= cos =(,...,m 
Beispiel 
Die Funktion sin 5 x in [-1, 1] soll durch ein Polynom vom Grad n=3 approximiert 


werden. Zur Verfügung stehen die Funktionswerte an den Nullstellen von Ts (x): 


x; | 951 588 0 .588 .951 
f; 997 -798 0 .798 .997- - 
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Maggetkarte einlesen (Block 1, 2) 
Programmbeginn 
Eingabe von: 


2 
2 
3 
4 
0 
1 
2 
3 
4 
5 
6 
7 
8 
9 
0 


1 


Ende der Koeffizienteneingabe 

Anzeige von: &9 0 
a, 1.1338968 
a2 0 
a3 -.1385336717 


Programm 4.5 Methode der kleinsten Quadrate 
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4.6 Der Algorithmus von Clenshaw 


Ist das Polynom 
pol (x) =a0To(x) +... +a,T„ (x) 


als Linearkombination von Tschebyscheff-Polynomen gegeben, so liefert das Programm den 
Wert des Polynoms an einer Stelle x = xo. 


Programminstruktionen 


Teen Jemen [mm ame 


1 | Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von ag, ..., An 


Ende der Koeffizienteneingabe 
Eingabe von x%9 


Ergebnisanzeige pol (xo) 
anderes Argument pol (xo) 


Eingabe von x; 
Ergebnisanzeige pol (x}) 


Registerinhalte 


Roo. -.., Rog: Programmzeiger 
Ror, Bei Rn+7: 90, +++, An 


Bemerkung 


Die Schritte 7 bis 9 der Programminstruktionen lassen sich beliebig oft wiederholen. 


Beispiel 
Das Polynom 
pol(x)=7 Tolx) +3 Tı(x) +6 Ta (x) +4 Ta(x) + Ta (x) 


soll an den Stellen xo = 1.5 und x} = - 1.5 ausgewertet werden. 
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Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von: a9 
a 
a 
a3 
a4 
Ende der Koeffizienteneingabe 


SW PRUD- DO - 


Eingabe von: xo 
Anzeige von: pol (xo} 
anderes Argument 
Eingabe von: x; 
Anzeige von: pol (x}) 
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4.7 De Castlejau 


Die Bernstein-Polynome vom Grad n 
B’ (A) = (7) 1-9; r=0,...,n 


bilden eine Basis für die Polynome bis zum Grad n im Intervall [O, 1]. Ein in Bernstein- 
Polynomen entwickeltes Polynom 


pol (A) = boBA(A) +... +b,BNA) 


heißt Bezier-Polynom, die Koeffizienten b; heißen Bezier-Punkte. Nach de Castlejau 
berechnet sich der Wert 


b,,... [A =), BEL 


i=r 


des Bezier-Polynoms vom Grad s-r zu den Bezier-Punkten b,, ..., b, an der Stelle A 
nach der Rekursion 


b,, s“ (1-X%) b,, s-ıtAbr41,...,s- 


“u 


Magnetkarte einlesen (Block 1) 
Programmbeginn 


Eingabe von r und s 


Eingabe von b,,... 


‚bs 


Ende der Koeffizienteneingabe 


Eingabe von A 


Ergebnisanzeige 


Registerinhalte 


Roo, ---, Roz: Programmzeiger 
Ros. Eee R,_r+8: b,, ae b, 
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Bemerkung 
Die b,, ..., b, werden „überschrieben’’; daher ist die Auswertung an nur einer Stelle A 
möglich. 
Beispiel 
Das Bezier-Polynom 
pol (A) =2BZ(A +3 BF(A) +4 BZ (A) +2 B3(M 


ist an der Stelle X = 1/3 auszuwerten. 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von: 


1 
1 
0 
0 
1 
2 
3 
4 
4 


Ende der Koeffizienteneingabe 
Eingabe von: X 
‚3333333333 


Anzeige von: bo, 1,2,3 2.888883889 


Programm 4.7 De Castlejau 


od FE LEL TS - 4a 42 RÜL 
ooi1 11 HA 1 1 In So Re u 
als 21 R=s 25 = ISO 42 STD 
ouE 43 STO a4 + - al 02 03 
co ana 42 STD 423 FÜL 
coS 42 5TO I u ol 01 
DI u u Pe Be +42 FÜL Ya = 
cr 21 F-& ci cd 43 FÜL 
203 42 =TO al 75 - jun] InIn] 
ca 01 0 = FÜL SI + 
all 2 DIN] wie] UT Fi 
cil 30 = 35 = 
DE 42 5STD 42 ST0O 
013 IE INS JS IS 
IB "& LEL E90 + 
15 au 1° 01 1 
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IE = 14 TI - 
42 =TD 44 42 FOL 
O4 og 05 aa il 
TE LEL ar SS + 
a4 [ek uB 1 = 
"a Rls 2 Il2 99 = 
a4 Er 102 42 370 
EIN u: 104 05 05 
43 FÜL 35 105 F3 Rl# 
DE CE 4: Is 05 05 
2 + [1 1c7 21 RES 


4.8 Bezier-Kurve 


Zur Approximation einer Funktion f(x) im Intervall [O, m] ist es günstig, das Intervall 

durch die Trennstellen k = 1,2, ....,m-1 in Segmente aufzuteilen und f(x} in jedem 
Segment durch ein Bezier-Polynom vom Grad n anzunähern. Dazu wird im Segment [k, k + 1) 
der Parameter A=x -k eingeführt; die Bezier-Punkte dieses Segments bezeichnet man mit 
Banks bak +7; ws Dak+n- 

Das Programm bestimmt den Wert der aus den Bezier-Polynomen zusammengesetzten 
Bezier-Kurve bez (x) an der Stelle x [O, m] als Wert des Bezier-Polynoms im Segment 

[k, k+1) an der Stelle A\=x-k nach dem Algorithmus von de Castlejau (siehe 4.7 


„De Castlejau’’). 
Baken 


6) k k+1 m 


Programminstruktionen 


Verfahren 


Magnetkarte einlesen (Block 1) 
Programmbeginn 
Eingabe von n und m sowie bo, ..-, Dnm 


Ende der Koeffizienteneingabe 


Eingabe von x 


Ergebnisanzeige 
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Registerinhalte 


Roo: ---, Roa: Programmzeiger 
Ro. Se Ram+1o: bo, PERF baum 


Bemerkung 


Die b; werden „überschrieben‘’; daher ist die Auswertung von bez (x) an nur einer Stelle x 
möglich. 


Beispiel 
Die im Intervall [O, 2] durch die Tabelle 


iloı23345%6 
bla ı 23413 


gegebene Bezier-Kurve vom Grad n =3 soll an der Stelle x = 1.5 ausgewertet werden. 


Magnetkarte einlesen (Block 1) 


Programmbeginn 
Eingabe von: n 


Eingabe von: 


1 
1 
3 
0 
1 
2 
3 
4 
5 
6 
7 
7 


Ende der Koeffizienteneingabe 
Eingabe von: x 
Anzeige von: bez (x) 
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Programm 4.8 Bezier-Kurve 
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4.9 Interpolation durch kubische Splines 


Gegeben seien n+1 Knoten 
(x. fx), i=0,...,n 


mit äquidistanten Stützstellen x;=xo tih, h fest. Stellt man zur Approximation der 
Funktion f an eine Bezier-Kurve s (siehe 4.8 „Bezier-Kurve‘’) die Forderung der zwei- 
maligen Differenzierbarkeit an den Trennstellen x; und gibt man zusätzlich die beiden 
Randbedingungen 


s’(xo) =s’(xn)=0, 
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so ist s dadurch eindeutig festgelegt. Eine solche Kurve heißt kubischer Interpolationsspline 
und mit diesen Randbedingungen natürlicher Spline. Das Programm berechnet die Koeffi- 
zienten der Polynome 


Px)satbiix-x)te; (x -x)?+d, (x -x,)? 
xEel[x;-1,x%:] i=1,...,n, 


aus denen sich s zusammensetzt. 


fo f, 


Anzeige 


Magnetkarte einlesen (Block 1, 2) 
Programmbeginn 


Eingabe der Nummer des ersten zu 
belegenden Speicherplatzes k >29 


Eingabe von f(xo}, ..., f(xn) 


Ende der Koeffizienteneingabe 
Eingabe von h,n und k 


Ergebnisanzeige 


Registerinhalte 


Roo, ---, Rog: Programmzeiger Re+2n+13 +, Rk+an: Pi, ---, En 

Re, ee Re+n: f(xo), „..7 f(xn) Re+3n+1: a Re+an: d;, 3 dn 
Resntiv--, Rk+2n! €1, :.., En Rek+4n+1r +, Rk+5n-ı: Zwischenergebnisse 
Bemerkung 


Das Programm berechnet kubische Interpolationssplines für bis zu n = 10 Segmente. 


Beispiel 
Die Funktion f(x) = sinnx soll im Intervall [O, 2] interpoliert werden. Es steht folgende 
Tabelle zur Verfügung 
x; 005 1 15 2 
tx) | 0 10-10 


Damit ist also n=4 und h=0.5 


Magnetkarten einlesen (Block 1, 2) 
Programmbeginn 
Eingabe von: k 
f(xo) 
fix, ) 
f(x2) 
f(x3) 
f(xa) 
Ende der Koeffizienteneingabe 


o 
PO mn POD-OnNDMD 


Eingabe von: h 


Anzeige von: 


0 
3 
0 
4 
1 
0 
6 
4 
0 


I 
I | 
-2ENG@ 


l 
20 
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Programm 4.9 Interpolation durch kubische Splines 
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5 Numerische Differentiation und Integration 


5.1 Numerische Differentiation 


Um näherungsweise die Ableitungen (X einer Funktion f an einer Stelle x zu bestimmen, 
liegt es nahe, f in der Umgebung von x durch ein Stützpolynom vom Grad n mit n+1 
äquidistanten Stützstellen x;=xo tih, i=0,...,n, anzunähern und dieses zu differen- 
zieren. Das Programm liefert Näherungen für die in der Tabelle aufgeführten Ableitungen 


9 = (ko +rh). 


Anzahl der Stützstellen Das Programm berechnet Näherungswerte für: 


fra 


fay2 


fa 


xXTXog 


Einzugeben sind außer den Stützwerten fo, ..., fa die Zahlen m:=n+1, r:= 
und k := Ordnung der Ableitung sowie h. 


Programminstruktionen 


Magnetkarte einlesen (Block 1, 2) 


Programmbeginn 
Eingabe von m, r,k, h, fo, ..-, fn 


Ende der Koeffizienteneingabe 


Ergebnisanzeige 
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Registerinhalte 


Roo, ---, Ros: Programmzeiger 


Bemerkung 


Werden die Zahlen m, r und k in einer Kombination eingegeben, die in der Tabelle nicht 
auftritt, so hält das Programm und der Rechner zeigt dies durch eine blinkende Anzeige an. 


Beispiel 


Gesucht ist eine Näherung für die 1. Ableitung der Funktion f(x) = sinx an der Stelle m/4. 
Gegeben ist die Tabelle 


x 0 n6 nS3 m 
sn, | 0 172 V32 1 
Es ist also 


ie n/4-0 
h n/6 


=15, k=1 und h=n/6. 


Magnetkarte einlesen (Block 1, 2) 


Programmbeginn 
Eingabe von: 


3.141592654 


‚5235987756 


1.732050808 


.8660254038 
7 
Ende der Koeffizienteneingabe 


Anzeige von: f,,, .7068616846 


Der exakte Wert ist .7071067812 
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Programm 5.1 Numerische Differentiation 
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5.2 Sehnentrapezsumme 
b 


Als Näherung für den Wert IK dx benutzt man die Sehnentrapezsumme Sn, die man 


a 
durch Aufteilen des Intervalls [a, b] in 2N Teilintervalle erhält. Es ist 


2 be 
(ftasswı+2 t(ari a) 
i=1 


Die Funktion f wird als Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 


b-a 
oN+1 


Sn” 


1. Die Funktionsvorschrift ist in Klammern einzuschließen. 
2. Für x ist RCL 00 zu setzen. 
3. Die Taste = darf nicht verwendet werden. 


4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


1 | Magnetkarte einlesen (Block 1) 

2 | Eingabe der Funktionsvorschrift f(x) 
3 | Programmbeginn 

4 | Eingabe von a,b und N 

5 | Ergebnisanzeige 


Registerinhalte 


Roo, :.., Ros: Programmzeiger 


Beispiel 


1 

Gesucht ist eine Näherung für I = | ax z mit N=3. 
1+x 
0 
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Magnetkarte einlesen (Block 1) 


Eingabe der Funktionsvorschrift 


Programmbeginn 
Eingabe von: 


Anzeige von: Sn .71847471236 


Es ist | = arctan 1 = .7853981634 


Programm 5.2 
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DE ÜE. 43 FEL 
44 SLUM On 05 
InIR jur] ES = 

AT DSE +43 RUL 
o4 ud le DE 
12 E SE = 

Te LEL 1 RrS 
12 C E LEL 


5.3 Romberg-Integration 


b 


Betrachtet man zur Approximation des Wertes ro dx die Folge der Sehnentrapez- 

a 
summen So, ..., Sn, 50 läßt sich die Konvergenz dieser Näherungsfolge durch wiederholte 
Extrapolation nach der Formel 


k=1,....N i=k-1,...,0 


beschleunigen. Das Programm liefert als Näherung für das Integral den Wert So, ...,n- 
Die Funktion f wird als Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift f{x) ist in Klammern einzuschließen. 
2. Für x ist RCL 00 zu setzen. 
3. Die Taste = darf nicht verwendet werden. 


4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Transen I Ten [rm ame 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift f(x) 
Programmbeginn 

Eingabe von a,b und N 
Ergebnisanzeige 
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Registerinhalte 


Roo, -.., Rog: Programmzeiger 

Rıo, sang Rn+1o0: So. PN, Sn 

Bemerkung 

Die So, ..., Sy werden durch die S; x „überschrieben”. 
Beispiel : 


Gesucht ist eine Näherung für I = | ex. mit N=3. 
1+x 


[e) 


Anmerkungen Eingabe 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift f{x) 


Programmbeginn 
Eingabe von: 


Anzeige von: $o,1,2,3 


Es ist | = arctan 1 = .7853981634 


.1853964459 
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Programm 5.3 Romberg-Integration 


uinie) gi FREI 102 Ss3 04 4 
cc [W ci 103 a4 43 FREI 
o02 EZ ü in) 104 sa 01 01 
cos D 1] 42 5TO 105 SE 52 & 
cc 3 oss 05 05 108 Sr’ 42 FOL 
nlnb] EZ ; z 107 Sa 01 01 
DIE 108 3 


Ti 


[Ina 
NINE} 
[nIwE-) 
oc 
cl 
Olz 
013 
ld 
015 


Ü1E 


109 
11C 


TI ıTı 


er 
ven 
7 
MIA 
ızı 
% 
De 
[8X 


mn 
ZT .- 
17 


Ir) 
Du vu 
NT 


ea TEN La 
Lu } 


[ee BEE EEE EEE SE EEE EEE RE GE EEG EEE EEG EG EFOEREFOSEEFG 


I CR IR 0 u 2 a Se 


O1” 42 =TD "a - 

"12 0 Ba di ci 1 

01a 4 -L 4 

ozn ü = as = 

Ü ri 5ER EST 
1) 24 [8 l2 03 
I BI + [1 1 

[ 43 RÜL Es 


Deu I 
EN ee ET Je dd 
1 
- 


Kt a a a pp pn pn pn 


De EN 


(Bu nu I au Wen We) 


I 


in 


f 
Tr 
a 
— 41 
FOAMNMTOCO 


I. 


2 
Dem] 

1) 
u 
peu} 


ı pn 


a1 IK us di 
= 4 


oe 1-4 


1 fa rd 


04542 FRI 
O4 1a 1u 
147 TE LEL 
4 26 
DT Sr 1 
a0 EETHW 


DRIN Barren 


Su 
n 


T 
DE Eu min 


—ı 
Er 
Tai 
4m“ 
Do Tu 


Lafer 


FE u 
ımını 
1 


140 5 Numerische Differentiation und Integration 


5.4 Das Eulersche Polygonzugverfahren 


Gegeben sei eine gewöhnliche Differentialgleichung erster Ordnung mit Anfangsbedingung 
vfß,ylıyla=y,. 


Gesucht ist y{b) = y„. Die Grundidee des Polygonzugverfahrens besteht darin, das 
Intervall [a, b] in n gleiche Teile zu teilen und die Lösungskurve y(x) durch den Strecken- 
zug mit den Ecken (x,,n,) zu ersetzen. Beim Eulerschen Polygonzugverfahren ist 


b-a 


Kk+17 x“r nn. X 7a 


b-a 
Mt ad, K=0,..,n=1. 


Die Funktion f wird als Unterprogramm eingegeben. 
Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift f(x, y) ist in Klammern einzuschließen. 

2. Für x ist RCL 00, für y RCL 01 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift f(x, y) 


Programmbeginn 


Eingabe von a,b, y, und n 


Ergebnisanzeige 
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Registerinhalte 


Roo, -.-, Rosa: Programmzeiger 


Beispiel 


Gegeben ist die Anfangswertaufgabe y’ = ge y(1)=2. In n= 10 Schritten soll eine 
Näherung für y(1.5) gefunden werden. 


Magnetkarte einlesen (Block 1) 


Eingabe der Funktionsvorschrift f(x, y) 


Programmbeginn 


Eingabe von: a 
b 


Anzeige von: 2.287646401 


Esist y=v/x?+3 und y(1.5) = 2.291287847 . 


‚Programm 5.4 | Das Eulersche Polygonzugverfahren 


Joo 76 LEL 012 42 STD 024 42 5STD N26 43 RÜL 
co 11 A a13 03 03 025 02 02 037° 02 02 
DlZ 91 R/S 014 43 REL D26 76 LEL sg 44 5LUM 
003 42 STD 015 04 04 0Z7 12 6 0133 00 00 
004 00 00 Me 75 - DEZE 71 SER 40 497 DSEZ 
0053 91 RS 017 43 RÜEL DEI 32 HE 04] 03 03 
ade 42 5STD olS 00 00 0306855 ud 12 68 

oor7 04 04 019 95 = G21 43 REL 043 43 FÜL 
ü0üs 91 R/S 02055  * ü32 0202 a44 01 01 
02 42 STD 021 43 RÜEL üS2 95 = 045 91 R3 
oo 01 01 022 03 03 üs4 44 SUM D4E 76 LEL 


coll 91 RFS 02395 = 0535 01 01 04733 x 
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5,5 Das Verfahren von Heun 


Gegeben sei eine gewöhnliche Differentialgleichung erster Ordnung mit Anfangsbedingung 
yafiy), v{a)=y,. 


Gesucht ist y({b) = y„. Das Programm bestimmt eine Näherung für y) Inn Schritten 
(für i=0,...,n-1) nach den Formeln 


nm + Fl, n) 
mit Fi lfothn) 
und fo=h-f(x,n), ı=h fs +h,n,+th°fo) 


mit h= >: und x;=a+tih. Die Funktion f wird als Unterprogramm eingegeben. Dabei 
ist folgendes zu beachten: 


1. Die Funktionsvorschrift f(x, y) ist in Klammern einzuschließen. 

2. Für x ist RCL 00, für y RCL 01 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Tieren Ten Frame an 


Magnetkarte einlesen (Block 1) 

Eingabe der Funktionsvorschrift f(x, y) 

Programmbeginn 

Eingabe von a, b, y, und n a 
b 
Ya 
n 

Ergebnisanzeige 
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Registerinhalte 


Roo: ---, Ros: Programmzeiger 


Beispiel 


Gegeben ist das Anfangswertproblem y’ = 5, y(1) =2. Inn = 10 Schritten soll eine 
Näherung für y(1.5) bestimmt werden. 


Anmerkungen Eingabe 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift f(x, y) 


Programmbeginn 
Eingabe von: a 


Anzeige von: 7,9 2.291263323 


Es ist y=/x?+3 und y(1.5) = 2.291287847. 
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| Programm 5.5 Das Verfahren von Heun 


r 


mE 3 le GE 


r 


TE LE I cl z 
ol Ai RA 35 = I45 42 570 25 + 
az 21 R-S5 SS + Ode 01 ul 42 RÜL 
[NIu RES 42 STD 42 FOL 4 42 FÜL l DS 
De nn] 0203 1b u Eu! Ei = 
AS 931 RS a5 = 4344 & 42 5STD 
[N TE Tr: 56 0 gen Sc ] 21 
In I E 04 051 E SE 
oa al FE NSE 2 
DuT b- Bae" See 12 USE 
BEN] 1 04 
c4l = 55 
OLE SE 
212 ur IT 7 £ 
old EIS OH sa E 
(15 43 FEL 453 i Ü: 
Ole c4 uP3 jT=AR] \ 1: 
or 25 + 0El 935 = +! 
D1E 42 FOL [u -Po 2 Ep In 
0193 os 05 E33 423 RÜL a 
Jzc ss = jRI- zT n 
url 4d= =TO EG 7055 + = 


5.6 Das klassische Runge-Kutta-Verfahren 


Gegeben sei eine gewöhnliche Differentialgleichung erster Ordnung mit Anfangsbedingung 
v-f,y), via=y,. 

Gesucht ist y(b) = y,. Das Programm liefert eine Näherung n, für y,, nach den Formeln 
Nm trFiku.n), i=0,..,n-1 

mit F= fo +2f,+2f,+f3) 


Pan = h fo 
und fo=hfix,n), fı=h‘f “truntz) 


2 
f2=hf xt Nr D f3=h-f(x,+th, n+f,). 


Dabei ist x;=ati'h, h> ae und n die Anzahl der Schritte. Die Funktion f wird als 


Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift f(x, y) ist in Klammern einzuschließen. 

2. Für x ist RCL 00, für y RCL 01 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 
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Programminstruktionen 


Team ______ [em TrmeT ame 


Magnetkarte einlesen (Block 1) 


Eingabe der Funktionsvorschrift f(x, y) 


Programmbeginn 
Eingabe von a, b, y, und n 


Ergebnisanzeige 


Registerinhalte 
Roo, ---, Roz: Programmzeiger 


Beispiel 


Gegeben ist das Anfangswertproblem y’ = 2 y(1)=2. Gesucht ist in n = 10 Schritten 
eine Näherung für y(1.5). 


Magnetkarte einlesen (Block 1) 
Eingabe der Funktionsvorschrift f(x, y) 
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A 1 


Programmbeginn 


Eingabe von: 


Anzeige von: n, 2.291287849 


Esist y= x?+3 und y(1.5) = 2.291287847 . 


Programm 5.6 Das klassische Runge-Kutta-Verfahren 


LER EYES 
Lv 


Zı Ja rc i 


4 
Ü 


u] 


RN Toın 


ıZı 
ma-.D 


Inc TE 55 JE 42 5STO il2 935 = 
ol Al 43 FÜL oe 07 07 103 44 5LM 
DZ 21 RS 04 04 url 39 + 104 00 00 
Oc2 42 5STD a5 = af 02 3 10542 RÜEL 
2o4 oo an 42 5TD ITZ 95 = Ile 07 07 
005 91 R-5 ur 72 44 ZUM 107 88 + 
CE 48 ES + ür’4 06 06 102 42 FÜL 
[Tu nF 06 rs 43 FELL 102 0505 
TE 35 = Te DT 07 En 
[UT IE B S 44 SLIM ürrT 55 +42 
OE DE ITS DE 
+2 FECL 79085 
a4 ud 20 42 REL 
Je] 05 = 
=) 
1 
1 
E 


O0 7 IE, 

43 FÜL E53 u 99 = 

ur 07 > 1 44 UM 

ES + 24 Er 1 11 

02: 204 i 43 RÜL 

I ge + al 95 = led 068 06 

Ei 43 FOL 231 42 370 125 42 5T0O 
ui 05 05 032 07 07 126 05 05 
43 Fü 35 = u32 55 = 127 Ar 152 
ca 03 42 5STD 04 0323 18: a3 02 
Eu Si 01 035 935 = 128 12 8 
42 3STD Fl SER ü3E 44 ZUM 120 43 RCL 
a4 04 33: ar 0 a 17 
FE LEL & x 135 423 RÜL “1 Res 
IE 8 43 RCL 2a ng F& LEL 
" ER 4 04 in 55 = EHE 
ER 2 5 = 101 02 58 
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5.7 Einschrittverfahren mit Schrittweitensteuerung 


Um Schwankungen des lokalen Diskretisierungsfehlers beim Rechnen mit konstanter Schritt- 
weite zu vermeiden, steuert das Programm bei der Lösung des Anfangswertproblems 
y’=flx,y}, yla) =y, die Schrittweite selbständig so, daß der lokale Diskretisierungsfehler 
annähernd konstant bleibt. Es benutzt dabei ein Paar von Runge-Kutta-Verfahren F,, Fo+ 1 
der Ordnung p bzw. p+1 und bestimmt die jeweilige Schrittweite nach der Formel 


€ 
Kama 
u N ho IFosı -Folte: 0.08P*' 


wobei € >O eine vorzugebende Toleranzschranke ist. Das Programm verwendet die beiden 
folgenden Verfahren F>, F>: 
p=2: nn t+Falx,n) 

mit Fass lforh) 


und fo=hfbo,n), h,hehfix,+th,n,+hf,) 


p=3: Nat als, N) 
mit F3= A + f} +4f,) 
und fo=h fix,n), fı,shfo+thhn,+thefo), 


h h 
ante, n+ 2 to+th)) : 


Dabei ist x; =x;_-1 + h,. Die Funktion f wird als Unterprogramm eingegeben. Dabei ist 
folgendes zu beachten: 


1. Die Funktionsvorschrift f(x, y) ist in Klammern einzuschließen. 

2. Für x ist RCL 00, für y RCL 01 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


2| Eingabe der Funktionsvorschrift f(x, y) 


3 | Programmbeginn 


4| Eingabe von a,b, y,, ho, €, P 


5 | Ergebnisanzeige 


Registerinhalte 


Roo, -.., Rog: Programmzeiger 
Rio, Rı1, Rıa: fo, fı, fa 


Bemerkung 


Soll ein anderes Paar von Runge-Kutta-Verfahren F,, Fa+1 verwendet werden, so sind diese 
als Unterprogramme D’ und D zu programmieren (die im Programmausdruck angegebenen 
Unterprogramme sind dann selbstverständlich wegzulassen). Dabei ist zur Auswertung der 
Funktion f der jeweilige x-Wert in Roo, der jeweilige y-Wert in Ro, zu speichern; h ist in 
Ros, X; in Ro2 und n, in Roa gespeichert. Die Funktionsvorschrift f(x, y) ist als Unter- 
programm x? zu programmieren. 


Beispiel 


Gegeben ist das Anfangswertproblem y’ = ” y(1)=2. Gesucht ist eine Näherung für y(1.5) 
mit ho = 0.1 und e = 10°°. 
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2 


Magnetckarte einlesen (Block 1, 2) 


Eingabe der Funktionsvorschrift f(x, y) 


Programmbeginn 


Eingabe von: a 


1.5 

2 

0.1 
0.00001 


Anzeige von: 2.291291758 


n 


Programm 5.7 Einschrittverfahren mit Schrittweitensteuerung 


Doc FE LEL 1 Rrs 42 43 RFÜEL JE 428 5STD 
Sol 1198 25 + a4 02 03 INT =T Bu u SE u rn 
oz 91 Res od 1 44 75 > Es FI 
2o2 42 5STD 35 = o4S 43 RCL les 71 SER 
cos ld 42 STO ode us 08 DEF 12 0" 
os 4 STo ur ur ü4? 35 = lDeE 95 = 
nu] UP WR 23 CR 43 42 5TO 062 So I=l 
or ag TE LBL 043 05 05 a 
DaIRb= 22 Ip osco Fe LEL a 43 REL 
ale) 3 FÜL [Si 22 LH ([r2 05 05 
cin us 0z SE 42 FEL ür2 95 = 
ci ES + 52 05 05 ür4d 25 + 
Liz 43 REL cS4 67 Eu IS 93% 
013 5 08 055 44 ZUM ‚ 00 0 
ul4 var USE ZE INWY " 088 
015 43 FÜL IST FT GE Ur 1-2 
de 02 03 055 44 SUM 17a 42 FEL 
ur 35 = 159 FE LEL su 0707 
c1lE z2 IH“ ED 24 CE [U Ich Base Sr E 
119 FF GE 0&1l 71 SER lSE 42 REEL 
NEO ZE LHE Oez 14 D 032 06 06 
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ge = Ss 05 124 

25 1 35 = 155 205 

55 42 STD len E» 

43 ou 00 187 r 

LIE, 42 57T 162 Fe: 

49 U2  uE 163 zu n 

45 43 FELL 1F0 ln . 

= a4 dd 171 211 85 x 

42 EI + 1rE el2 04 4 

Din 423 FOL >23 2i3 85 + 

22 1. 15 05 74 214 43 ROLL 

45 EI“ 75 sis 1i il 

E = 42 FÜL TE ZiE 85 + 

& 25 [PT Tr ziT 42 FEL 

3 n 5 = TE „ eis 10 10 

ae 42 ST 79 | 213 54 5 

& = Oi 0 53] :ER 220 55 = 

43 RCL 42 5 = fr i £ 
il 


L 
ic 
m in 
ee ER 


NOSWrEOrD 


027 
ll 
oO 
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Fa rc 


TE EN PLN LS PLN N en 


eo 


RE er 
+ 


Pe 

Fa 

a 
u Wi 


EN Fe 


Pet 
De 
De 
De} 
2 
De} 


r“ - ide 05 u ia ı0 D* 

42 ROL 147 &1 GTD 54 il 
EN) JS 142 Z2 IMY ‘ Ü 
12 958 = 142 FE LEL FOL 
1lü 7 GE 150 44 SuM 10 
1il 25 CLRF 121 423 FÜEL SS + 
liz 4 122 04 04 42 FÜL 
112 u 153 21 ES 11 ıi 
il4 4 154 TE LEL 4 
115 0 155 14 I s5 = 
llE 58 13652 & DEE 
117 2 12742 RÜL 4 
ilS 7 155 5 32 RTH 
113 8 1234 " LEL 
lzu 0 4 180 BE 
121 u 1&1 
1 l&z 
1E 18: 


5.8 Die Mittelpunktsregel 


Die Mittelpunktsregel ist ein Mehrschrittverfahren zur Lösung der gewöhnlichen 
Differentialgleichung erster Ordnung 


yaf,y), 


bei dem vor Beginn der Rechnung zwei Startwerte y(xo) = yo und y(xı) =y, bekannt 
sein müssen. Eine Näherung n für das gesuchte y(b) = y, bestimmt sich dann nach den 
Formeln 


Pa nu 2 ES A TPerZ Per EZ 
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Dabei ist x;=a + ih. Die Länge des Intervalls [a, b] muß ein ganzzahliges Vielfaches der 
Schrittweite h = x} - xo sein, also 


b-a=n-h mit nen. 
Die Funktion f wird als Unterprogramm eingegeben. Dabei ist folgendes zu beachten: 


1. Die Funktionsvorschrift f(x, y) ist in Klammern einzuschließen. 

2. Für x ist RCL 00, für y RCL 01 zu setzen. 

3. Die Taste = darf nicht verwendet werden. 

4. Die Eingabe der Funktionsvorschrift ist mit INV SBR abzuschließen. 


Programminstruktionen 


Tran Transen [ram rn 


Magnetkarte einlesen (Block 1) 

Eingabe der Funktionsvorschrift f(x, y) 

Programmbeginn 

Eingabe von a, b, yo, Yı. h a 
b 
Yo 
Yı 
h 

Ergebnisanzeige 


Registerinhalte 


Roo; ---, Ros: Programmzeiger 


Beispiel 


Gegeben ist die Differentialgleichung y’ = 5 mit den Startwerten y(1)=2 und 
y(1.05) = 2.025. Gesucht ist eine Näherung für y(1.5); es ist h = 1.05 - 1 = 0.05. 
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Magnetkarte einlesen (Block 1) 


Eingabe der Funktionsvorschrift f(x, y) 


Programmbeginn 


Eingabe von: a 
b 


1.5 
2 
2.025 


Anzeige von: 2.291400512 ) 


Esist y=x?+3 und y(1.5) = 2.291287847 . 


Programm 5.8 Die Mittelpunktsregel 


ocl Te LEL +43 FÜL 
ol 11 A u 03 
clz 921 R#S = 
clz 32 STO 42 ERl 
004 00 00 co 0 
cos 931 R#S 42 STD 
ole 42 5TO O3 02 
Du De a 1 Fo ne ar IS 
co a1 RS 05 05 
03 42 STD l2E EB 
ol0 02 03 43 FRÜL 
011 91 R-s ci 01 
Ol2 42 STD 1 RrS 
12 ci 0 r& LEL 
Ba ME 


old 91 Res 
015 42 570 
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Verzeichnis der behandelten Probleme 


DC BEL EEE 


Anfangswertaufgaben erster Ordnung 


Approximation 


Auswerten von Polynomen 


Differentiation 
— von Polynomen 


Eigenvektoren 


Eigenwerte 


Fixpunkte 


Integration 


Interpolation 
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Problem Programm 


Lineare Gleichungssysteme 


— symmetrische 


— überbestimmte 


Nullstellen 


— von Polynomen 


Taschenrechner-Literatur 


Aus der Reihe: Anwendung programmierbarer Taschenrechner 

Die Reihe „Anwendung programmierbarer Taschenrechner’’ bietet den Benutzern 
dieser Rechner eine breite Palette von Aufgabenstellungen aus den Anwendungs- 
bereichen der Natur- und Wirtschaftswissenschaften an, für die Rechnerprogramme 
zur numerischen Lösung dargestellt werden. 


Helmut Alt 

Angewandte Mathematik — Finanzmathematik — Statistik — Informatik für UPN- 
Rechner. 1979. VIIl, 162 S. DIN C 5 (Anwendung programmierbarer Taschenrech- 
ner, Bd. 1). Kart. 


Peter Kahlig 

Mathematische Routinen der Physik, Chemie und Technik für AOS-Rechner — Teil 1. 
Mit 71. Abb., 129 Beispielen und 34 Tabellen. 1979. VI, 178S. DINC 5 (Anwen- 
dung programmierbarer Taschenrechner, Bd. 3/1). Kart. 


Harald Nahrstedt 

Statik — Kinematik - Kinetik für AOS-Rechner. Mit 30 vollständigen Programmen 
und 140 Abb. 1980. ca. 120 S. DIN C 5 (Anwendung programmierbarer Taschen- 
rechner, Bd. 4). Kart. 


Aus der Reihe: Programmieren von Taschenrechnern 

Mit diesen Büchern werden dem im Programmieren unerfahrenen Leser Kenntnisse 
über den Umgang mit programmierbaren Taschenrechnern vermittelt. Jeder Band ist 
auf bestimmte Rechnertypen zugeschnitten. 


Paul Thießen 

Lehr- und Übungsbuch für die Rechner HP-29C/HP-19C und HP-67/HP-97. Hrsg. 
von Hans H. Gloistehn. 1980. VIII, 153 S. 12 X 19,5 cm (Programmieren von 
Taschenrechnern, Bd. 4). Kart. 


Hans-Joachim Ludwig 
Programmoptimierung für Taschenrechner (AOS). 1979. X, 102S. 12X 19,5 cm 
(Programmieren von Taschenrechnern, Bd. 5). Kart. 


Info-Gutschein 
Bitte informieren Sie mich (uns) ständig über ihre 
Neuerscheinungen auf dem Gebiet: 
DO Taschenrechner OD Mikrocomputer 
Ich (wir) besitze(n) folgendes Gerät: 


TR: 


uC: 


Hauptanwendungsgebiete des TR bzw. uC: 


Diese Karte entnahm (en) ich (wir) dem Buch: 
Kahmann, Anwendung programmierbarer TR, Bd.5 


Meine (unsere) Buchhandlung: 


Gleichzeitig bestelle(n) ich (wir) folgende Bücher: 


Autor und Titel 


Schumny, TR + „C-Jahrbuch 1981 


Böhm /Gose, Numerische Mathematik 


Anschrift: 


Beruf/Branche 


Datum Unterschrift 


Lieber Leser! 


Wenn Sie Interesse haben, aktiv an der 
Weiterentwicklung unseres Literaturprogramms 
zum Bereich TR + uC mitzuarbeiten, z.B. durch 
Veröffentlichung ausgetesteter Programme zu 
bestimmten Anwendungsgebieten, dann schreiben 
Sie uns. 


Wir freuen uns über Ihre Nachricht und werden 
uns umgehend mit Ihnen in Verbindung setzen. 


Mit freundlichen Gruß 
Lektorat Fachbuch 
gez. Niclas 


Bitte 
mit 


50 Pf. 
freimachen 


Antwort 


Friedr. Vieweg & Sohn 
Verlagsgesellschaft mbH 


Postfach 5829 


D-6200 Wiesbaden 1 


Anwendung programmierbarer Taschenrechner 


Diese Reihe bietet den Benutzern programmierbarer Taschenrechner eine reich- 
haltige Palette von Aufgabenstellungen aus den Anwendungsgebieten der Natur- 
und Wirtschaftswissenschaften an, für die Programme zur numerischen Lösung 
entwickelt werden. 


Jeder Band behandelt ein in sich abgeschlossenes Themengebiet: Nach einer kurzen 
Einführung in die Theorie der jeweiligen Problemstellung wird der Lösungsalgorith- 
mus entwickelt, das Programm dargestellt und kommentiert. 


Neben der direkten Nutzung der hier veröffentlichten Programme unterstützt diese 
Reihe den Leser wirkungsvoll bei der Ausarbeitung eigener Programmvarianten. 


Band 5: Numerische Mathematik 
Programme für den TI 59 


von Jürgen Kahmann 


Band 5 enthält 44 Programme zur numerischen Mathematik: Matrizen / Lineare 
Gleichungen und Ungleichungen / Iteration / Interpolation und diskrete Approxi- 
mation / Numerische Differentiation und Integration. 


Jürgen Kahmann ist Mitarbeiter am Institut für Angewandte Mathematik 
der Technischen Universität Braunschweig. 


ISBN 3-528-04171-4 


